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An extensive set of field measurements have been collected throughout the continental margin of the northeastern 
U.S. from 2004 to 201 1 to develop and validate ocean color satellite algorithms for the retrieval of the absorption 
coefficient of chromophoric dissolved organic matter ( Ocdom ) and CDOM spectral slopes for the 275:295 nm and 
300:600 nm spectral range ( S 27 5 :2 95 and S 30 o:6oo )• Remote sensing reflectance (R re ) measurements computed from 
in-water radiometry profiles along with Qcdom(^) data are applied to develop several types of algorithms for the 
SeaWiFS and MODIS-Aqua ocean color satellite sensors, which involve least squares linear regression of OcdomOs) 
with ( 1 ) R rs band ratios, (2) quasi-analytical algorithm-based (QAA-based) products of total absorption coefficients, 
(3) multiple R„ bands within a multiple linear regression (MLR) analysis, and (4) diffuse attenuation coefficient 
(IQ). The relative error (mean absolute percent difference; MAPD) for the MLR retrievals of Ocdom( 275), 
Ocdom( 355), a C DOM(380), o C dom(4 12) and cicdom(443) for our study region range from 20.4 to 23.9% for MODIS- 
Aqua and 27.3-30% for SeaWiFS. Because of the narrower range of CDOM spectral slope values, the MAPD for the 
MLR S 2 75 : 295 and S 30 o :6 oo algorithms are much lower ranging from 9.9% and 9.1% for SeaWiFS, respectively, and 
8.7% and 9.7% for MODIS, respectively. Multi-year, seasonal and spatial MODIS-Aqua and SeaWiFS distributions of 
Ocdom, S 275295 and S 300 : 6 oo processed with these algorithms are consistent with field measurements and the process- 
es that impact CDOM levels along the continental shelf of the northeastern U.S. Several satellite data processing fac- 
tors correlate with higher uncertainty in satellite retrievals of o C dom, $ 275.295 and S 30 o:6oo within the coastal ocean, 
including solar zenith angle, sensor viewing angle, and atmospheric products applied for atmospheric corrections. 
Algorithms that include ultraviolet R re bands provide a better fit to field measurements than algorithms without the 
ultraviolet bands. This suggests that satellite sensors with ultraviolet capability could provide better retrievals of 
CDOM. Because of the strong correlations between CDOM parameters and DOM constituents in the coastal ocean, 
satellite observations of CDOM parameters can be applied to study the distributions, sources and sinks of DOM, 
which are relevant for understanding the carbon cycle, modeling the Earth system, and to discern how the Earth 
is changing. 

Published by Elsevier Inc. 


1. Introduction 

Chromophoric dissolved organic matter (CDOM) represents the op- 
tically active fraction of DOM in natural waters. CDOM absorption is 
characterized as an exponential decrease in absorption from ultraviolet 
(UV) to visible wavelengths. CDOM absorption coefficients, Ocdom(A). 
and the spectral slope coefficient (S), a parameter that quantifies the 
exponential absorption decrease with increasing wavelength, have 
been shown to vary with type and source of CDOM (Blough & Del 
Vecchio, 2002; Bricaud, Morel, & Prieur, 1981 ). Terrestrial plant matter 
is considered to be the primary source of CDOM to the coastal ocean 
(Del Castillo, Coble, Morell, Lopez, & Corredor, 1999; Del Vecchio & 
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Blough, 2004). For example, Hernes and Benner (2003) found a strong 
correlation between dissolved lignin phenols (compounds derived 
from vascular plants) and a CDO M(350) within the Mississippi River 
plume. Such studies demonstrate that at least within the coastal ocean 
Ocdom may be useful as a tracer of terrigenous DOM. However, biological 
processes such as phytoplankton growth, zooplankton grazing and mi- 
crobial activity can also contribute marine-derived CDOM to continental 
margins and pelagic ocean (Andrew, Del Vecchio, Subramaniam, & 
Blough, 2013; Nelson & Siegel, 2002; Rochelle-Newall & Fisher, 2002; 
Steinberg, Nelson, Carlson, & Prusak, 2004). Characteristics of CDOM 
absorbance can be applied to trace the mixing of various water masses 
(e.g., Stedmon, Osburn, & Gragh, 2010). Physical processes that promote 
vertical mixing such as winter seasonal mixing, upwelling, and storms 
can introduce CDOM to the surface ocean. Deep convective mixing of 
the ocean by intense tropical storms such as hurricanes can also transport 
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CDOM from depth to the surface (Hoge & Lyon, 2002). Furthermore, re- 
cent work has shown that CDOM is strongly correlated to the apparent 
oxygen utilization and therefore could potentially be used as a tracer of 
biogeochemical processes and ocean circulation in the global ocean 
(Nelson, Siegel, Carlson, & Swan, 2010; Swan, Siegel, Nelson, Carlson, & 
Nasir, 2009). Nelson and Siegel (2013) recently reviewed the distribution 
and dynamics of CDOM in the global ocean. CDOM can dominate the 
inherent light absorption at blue wavelengths in surface waters of the 
coastal (20-70% at 440 nm; Del Vecchio & Subramaniam, 2004; Pan, 
Mannino, Russ, & Hooker, 2008) and pelagic ocean (>50% at 440 nm; 
Siegel, Maritorena, Nelson, Hansell, & Lorenzi-Kayser, 2002). Because 
CDOM is ubiquitous and a dominant light absorbing constituent in the 
ocean, accurate retrieval of CDOM absorption is a prerequisite for apply- 
ing ocean color remote sensing data to quantity other optically active 
ocean constituents such as chlorophyll a or rate processes such as phyto- 
plankton productivity and optical properties. The impact of CDOM 


absorption on photosynthetically active radiation (PAR) as well as the 
spectral quality of sunlight within the ocean's euphotic zone can influ- 
ence primary production, primarily in regions where phytoplankton are 
limited by light such as in higher latitudes or in turbid waters. Seasonal 
vertical stratification isolates DOM at depth from sunlight, which can 
degrade CDOM, resulting in DOM with greater chromophoric content at 
depth than at the surface (Nelson & Siegel, 2002; Vodacek, Blough, 
DeGrandpre, Peltzer, & Nelson, 1997). The loss of CDOM from surface wa- 
ters can be related to the photochemical production of dissolved inorganic 
carbon (e.g., Johannessen, Miller, & Cullen, 2003; Miller & Zepp, 1995). 

Dissolved organic carbon (DOC) in the ocean constitutes one of the 
largest pools of organic carbon in the biosphere (700 x 10 15 g C; Field, 
Sarmiento, & Hales, 2007) and nearly all (>97%) of the organic carbon 
in the ocean (Benner, 2002). Over time, a change in the balance between 
production and remineralization of DOC may impact the amount of C0 2 
in the ocean and ultimately the atmosphere. Hypothetically, an increase 
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Fig. 1. Map of the study region and station locations sampled for (a) algorithm development, (b) validation within the coastal regions of the northwestern Atlantic Ocean, which include 
the Gulf of Maine (GoM), Georges Bank (GB), and Middle Atlantic Bight (MAB), and (c) time series and monthly composite comparisons with the satellite data. The symbol colors in (a) and 
(b) represent the season sampled: blue = winter, green = spring, red = summer, orange = fall. The symbol shapes in (b) represent the satellite sensor applied in the algorithm valida- 
tion analysis: circle = MODIS-Aqua and triangle = SeaWiFS. The gray lines on (a) and (b) maps represent the approximate bathymetry contours of 20, 60, 100 and 500 m. The study re- 
gions are further subdivided into sub-regions based on bathymetry as indicated in the map legend and latitude (39° N) to distinguish the northern (NMAB) and southern MAB (SMAB) 
areas. The inset map shows the station locations from the lower Chesapeake Bay, mouth, plume and shelf waters for the time series comparison of field measurements and satellite 
retrievals of CDOM properties. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 
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in the oxidation of the DOC pool in the ocean (to C0 2 ) of 1% annually 
would exceed the amount of C0 2 released from fossil fuel combustion 
(Hedges, 2002). Higher rates of photochemical degradation of CDOM 
(e.g., from increasing UV radiation to the ocean surface) and microbial 
respiration (related to increasing ocean temperatures) are processes 
that can increase the net oxidation of DOC in the ocean. Despite com- 
prising a relatively small fraction of DOM, CDOM strongly correlates 
with DOC in rivers, river plumes, estuaries and continental margins 
(e.g., Del Castillo et al„ 1999; Ferrari, Dowell, Grossi, & Targa, 1996; 
Mannino, Russ, & Hooker, 2008; Spencer, Butler, & Aiken, 2012; 
Stedmon, Markager, & Kaas, 2000; Vodacek et al., 1997), but not in 
pelagic ocean waters (Nelson & Siegel, 2002). The strong relationships 
between CDOM absorption and DOM in coastal waters provide a means 
to quantify DOC as well as sources and sinks of DOM from satellite re- 
trievals of CDOM absorption. 

Numerous bio-optical algorithms have been developed to retrieve 
CDOM absorption or closely related products from ocean color satellite 
observations including SeaWiFS (Sea-viewing Wide Field-of-view 


Sensor), MODIS (MODerate resolution Imaging Spectroradiometer) on 
the Aqua satellite, and MER1S (MEdium-spectral Resolution Imaging 
Spectrometer). Many of these inversion models estimate the absorption 
coefficient of CDOM and detrital (non-pigmented) particles as a single 
parameter (a CDM ), because CDOM and detritus have similar spectral re- 
sponses in the visible spectrum (Bricaud, Ciotti, & Gentili, 2012; Carder, 
Chen, Lee, Hawes, & Kamykowski, 1999; Doerffer & Schiller, 2007; Hoge, 
Wright, Lyon, Swift, & Yungel, 2001 ; Lee, Arnone, Hu, Werdell, & Lubac, 
2010; Lee, Carder, & Arnone, 2002; Maritorena, Siegel, & Peterson, 2002; 
Siegel, Maritorena, Nelson, & Behrenfeld, 2005; Siegel, Maritorena, 
Nelson, Behrenfeld, & McClain, 2005; Tilstone et al., 2012; Werdell 
et al., 2013). In oligotrophic ocean waters such as the Sargasso Sea, 
non-pigmented particles account for only -12% of CDM absorption 
(Siegel et al., 2002). Recent work in coastal ocean waters indicates 
that non-pigmented particles account for <20% of total light absorption 
at 443 nm along the U.S. Mid-Atlantic coast (Pan et al., 2008). The vari- 
ous inversion models apply remote sensing reflectance ( R „) for several 
wavelength bands to simultaneously derive multiple parameters 


Table 1 

List of sampling cruises and measurements collected for the naperian absorption coefficient of colored dissolved organic matter (a CD0M ) algorithm development or validation. 


Region 

Cruise 

Dates 

Measurements 

Stations 3 

CDOM 

AOPs 

Lower CB b 

CBM01 

5 May 2004 

CDOM c 

2, 0,0 


Lower CB 

CBM02 

5 July 2004 

CDOM 

4, 0,0 


Lower CB 

CBM03 

1 September 2004 

CDOM 

4, 0,1 


Lower CB 

CBM04 

15 October 2004 

CDOM 

4, 0,0 


Lower CB 

CBM05 

1 5 November 2004 

CDOM 

4, 0,2 


Lower CB 

CBM06 

10 January 2005 

CDOM 

4, 0,0 


Lower CB 

CBM07 

26 May 2005 

CDOM 

3,0,1 


Lower CB 

CBM08 

21 June 2005 

CDOM 

3,0,0 


Lower CB 

CBM09 

19 August 2005 

CDOM 

4, 0,0 


Lower CB 

CBM 10 

23 September 2005 

CDOM 

4, 0,2 


Lower CB 

CBM 11 

30 January 2006 

CDOM 

3,0,0 


Lower CB 

CBM12 

12 April 2006 

CDOM 

4, 0,0 


CB Mouth & Plume 

D01 

27 May 2005 

CDOM 8i AOPs d 

4, 2,1 

4,2, 1 

CB Mouth & Plume 

D02 

3 November 2005 

CDOM & AOPs 

6, 6,0 

6, 6,5 

CB Mouth & Plume 

D03 

6 September 2006 

CDOM & AOPs 

5,2,0 

2,2,0 

CB Mouth & Plume 

D04 

28 November 2006 

CDOM & AOPs 

6,3,0 

3,3,0 

CB Mouth & Plume 

PL6 

19 March 2007 

CDOM 

4, 0,2 


CB Mouth & Plume 

COI1 

23 April 2007 

CDOM 

4, 0,2 


CB Mouth & Plume 

COI2 

3 July 2007 

CDOM 

4, 0,2 


CB Mouth & Plume 

COB 

16 August 2007 

CDOM 

4, 0,0 


Southern MAB e 

B01 

30 March-1 April 2005 

CDOM & AOPs 

14, 5, 0 

5, 5,0 

Southern MAB 

B02 

26-30 July 2005 

CDOM & AOPs 

29, 18,0 

18, 18,4 

Southern MAB 

B03 

9-12 May 2006 

CDOM 

26, 0, 4 


Southern MAB 

B04 

2-6 July 2006 

CDOM & AOPs 

31,30,0 

30, 30, 3 

Southern MAB 

B05 

31 Oct.-2 Nov. 2006 

CDOM 

16, 0, 1 


GoM f 

BIOD01 

26-30 April 2007 

CDOM & AOPs 

14, 9, 0 

9, 9,0 

GoM 

BIOD02 

26-28 May 2007 

CDOM & AOPs 

17, 14,0 

14, 14,3 

GoM & MVCO 8 

BIOD03 

6-8 June 2007 

CDOM & AOPs 

13, 11,0 

11, 11,0 

NYB 1 ' 

OCV1 

5-9 May 2007 

CDOM & AOPs 

20, 20, 0 

20, 20, 8 

NYB 

OCV2 

10-14 Nov. 2007 

CDOM & AOPs 

22, 16,0 

16, 16,4 

NYB 

OCV3 

21-23 July 2008 

CDOM & AOPs 

8, 8,0 

8, 8,0 

NYB 

OCV5 

19-21 May 2009 

CDOM & AOPs 

11, 11,0 

11, 11,3 

MAB, GoM, GB 1 

CV1 

17-28 August 2009 

CDOM 

51,0, 10 


MAB, GoM, GB 

CV2 

3-20 Nov. 2009 

CDOM 

69, 0, 19 


MAB, GoM, GB 

CV3 

1-17 February 2010 

CDOM 

45, 0, 4 


MAB, GoM, GB 

CV4 

26 May-9 June 2010 

CDOM 

79, 0, 14 


MAB, GoM, GB 

CV5 

6-21 Nov. 2010 

CDOM 

69, 0, 10 


MAB, GoM, GB 

CV6 

2-21 June 2011 

CDOM 

93, 0, 9 



a Number of stations measured or applied for CDOM algorithm development and validation (n, n, n), respectively and number of stations where AOP profiles were collected, used for 
CDOM algorithm development, and for validation of satellite remote sensing reflectance (n, n, n). For validation, only stations meeting the validation criteria for either MODIS-Aqua or 
SeaWiFS are listed. Blank cells indicate that no measurements were collected or applied. For a sub-set of stations, a second set of CDOM samples and series of AOP profiles were collected 
or conducted within approximately one hour of the first set of measurements. 
b Lower CB — Old Dominion University's Chesapeake Bay mouth hydrography transect north of the Chesapeake Bay bridge and tunnel system. 
c CDOM — spectral absorption coefficients of colored dissolved organic matter. 

d AOPs — apparent optical properties including spectral water-leaving radiances, remote sensing reflectances and light attenuation coefficients. 
e Southern MAB — Delaware Bay and Chesapeake Bay mouths, plume and adjacent continental margin of the southern Middle Atlantic Bight (MAB). 
f CoM — Gulf of Maine. 

8 MVCO — area surrounding the Martha's Vineyard Coastal Observatory (MVCO) south of Cape Cod, Massachusetts. 
h NYB — lower Hudson River, Hudson-Raritan Bay and plume and adjacent continental margin of the New York Bight (NYB). 

1 GB — Georges Bank. 
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including a CDM , a ph (phytoplankton absorption coefficient) or chloro- 
phyll a (Chi), and the particulate backscatter coefficient. Satellite prod- 
ucts are also derived through empirical methods that relate a CDOM to R rs 
band-ratios of various ocean color bands (e.g., Kahru & Mitchell, 2001; 
Shanmugam, 2011). Various R rs band-ratio algorithms [K r s(^i)/Krs(^ 2 )] 
have been developed to retrieve o C dom within coastal regions such 
as R rs (443 nm)/R re (555 nm) (D'Sa & Miller, 2003; D'Sa, Miller, & Del 
Castillo, 2006; Del Castillo & Miller, 2008; Johannessen et al., 2003; 
Kahru & Mitchell, 2001; Mannino et al., 2008; Naik, D'Sa, Grippo, 
Condrey, & Fleeger, 2011) and global ocean (Shanmugam, 2011; Tiwari 
& Shanmugam, 2011). Retrieval of a CDM using inversion algorithms yields 
reasonable results in open ocean regions (Siegel, Maritorena, Nelson, 
Behrenfeld, & McClain, 2005; Siegel et al., 2002). However, such algo- 
rithms typically do not work well in coastal waters due to the optical 
complexity of turbid coastal waters (high levels of CDOM, colored detrital 


particles, and phytoplankton; e.g., Aurin & Dierssen, 2012) as well as sim- 
plified parameterization (e.g., single S to represent marine waters global- 
ly; Maritorena et al., 2002; Magnuson, Harding, Mallonee, & Adolf, 2004; 
Morel & Gentili, 2009). Recent approaches have attempted to separate 
a CDOM from Qcdm (Dong, Shang, & Lee, 2013; Matsuoka, Hooker, Bricaud, 
Gentili, & Babin, 2013; Tilstone et al., 2012) and to derive the CDOM spec- 
tral slope in various spectral regions (Fichot et al., 2013; Mannino et al., 
2008). Despite the many algorithms developed to retrieve Ocdom and 
a CD M from satellite data only a few have undergone rigorous validation 
involving direct comparison of field measurements with coincident satel- 
lite data. The primary limitation to rigorous validation is the lack of suffi- 
cient data of coincident field measurements and satellite observations 
that are independent from the data used to develop the algorithm. 

Accurate retrievals of cicdom(A) are needed to improve satellite esti- 
mates of Chi. The intense and spectrally coincident absorption by CDOM 


Table 2 

Spatial and temporal range of in situ CDOM measurements. Mean +/— sd. 


Sub-region 

Season 

Qcdom3 55 


Qcdo/vj412 


5 275:295 


S 300:600 


n 



Mean 

±sd 

Mean 

±sd 

Mean 

±sd 

Mean 

±sd 


Gulf of Maine 

<50 ill 

Spring 

1.033 

0.750 

0.393 

0.295 

0.0213 

0.0025 

0.0180 

0.0008 

6 


Summer 

1.470 

0.202 

0.574 

0.083 

0.0182 

0.0009 

0.0169 

0.0004 

11 

50-100 m 

Winter 

1.437 

0.554 

0.562 

0.221 

0.0185 

0.0016 

0.0170 

0.0005 

14 


Spring 

0.419 

0.049 

0.171 

0.023 

0.0236 

0.0010 

0.0176 

0.0010 

8 


Summer 

0.563 

0.155 

0.226 

0.063 

0.0220 

0.0016 

0.0175 

0.0005 

6 


Fall 

0.632 

0.178 

0.240 

0.064 

0.0233 

0.0017 

0.0182 

0.0011 

18 

100-1000 m 

Winter 

0.571 

0.376 

0.219 

0.142 

0.0237 

0.0034 

0.0185 

0.0008 

6 


Spring 

0.332 

0.080 

0.131 

0.034 

0.0253 

0.0014 

0.0186 

0.0010 

12 


Summer 

0.302 

0.045 

0.121 

0.018 

0.0264 

0.0009 

0.0191 

0.0004 

25 


Fall 

0.384 

0.107 

0.146 

0.041 

0.0264 

0.0016 

0.0192 

0.0010 

75 

Georges Bank 

<50 in 

Winter 

0.261 

0.019 

0.109 

0.007 

0.0261 

0.0005 

0.0189 

0.0003 

5 


Summer 

0.369 

0.065 

0.159 

0.033 

0.0256 

0.0016 

0.0179 

0.0009 

14 


Fall 

0.326 

0.026 

0.136 

0.018 

0.0256 

0.0008 

0.0187 

0.0006 

3 

50-75 m 

Winter 

0.279 

0.021 

0.118 

0.007 

0.0259 

0.0002 

0.0185 

0.0004 

4 


Summer 

0.304 

0.040 

0.122 

0.020 

0.0275 

0.0011 

0.0191 

0.0008 

28 


Fall 

0.319 

0.038 

0.132 

0.018 

0.0261 

0.0006 

0.0186 

0.0006 

6 

75-200 m 

Winter 

0.265 

0.061 

0.105 

0.021 

0.0266 

0.0012 

0.0189 

0.0007 

10 


Summer 

0.283 

0.057 

0.110 

0.026 

0.0283 

0.0018 

0.0198 

0.001 1 

23 


Fall 

0.240 

0.099 

0.098 

0.040 

0.0296 

0.0047 

0.0194 

0.0008 

16 

Northern MAB 

<25 m 

Winter 

0.477 

0.103 

0.199 

0.036 

0.0235 

0.0012 

0.0176 

0.0004 

4 


Spring 

1.099 

0.410 

0.429 

0.164 

0.0206 

0.0028 

0.0180 

0.001 1 

14 


Summer 

0.650 

0.197 

0.260 

0.088 

0.0226 

0.0028 

0.0183 

0.0008 

20 


Fall 

0.537 

0.164 

0.218 

0.058 

0.0242 

0.0018 

0.0182 

0.0004 

15 

25-60 m 

Winter 

0.358 

0.081 

0.151 

0.032 

0.0253 

0.0015 

0.0180 

0.0005 

14 


Spring 

0.557 

0.263 

0.214 

0.105 

0.0245 

0.0029 

0.0194 

0.0013 

24 


Summer 

0.489 

0.170 

0.186 

0.067 

0.0256 

0.0023 

0.0192 

0.0012 

35 


Fall 

0.405 

0.081 

0.167 

0.032 

0.0255 

0.0011 

0.0183 

0.0004 

35 

60-100 m 

Winter 

0.239 

0.068 

0.105 

0.031 

0.0276 

0.0019 

0.0183 

0.0007 

5 


Spring 

0.284 

0.034 

0.103 

0.015 

0.0284 

0.0009 

0.0211 

0.0007 

4 


Summer 

0.285 

0.088 

0.106 

0.036 

0.0292 

0.0022 

0.0207 

0.0015 

23 


Fall 

0.312 

0.052 

0.129 

0.024 

0.0266 

0.0013 

0.0186 

0.0006 

15 

100-1000 m 

Winter 

0.190 

0.083 

0.080 

0.034 

0.0291 

0.0036 

0.0192 

0.0014 

3 


Summer 

0.203 

0.033 

0.074 

0.015 

0.0315 

0.0020 

0.0209 

0.001 1 

12 


Fall 

0.247 

0.052 

0.103 

0.023 

0.0282 

0.0012 

0.0188 

0.0005 

6 

Southern MAB 

<25 m 

Winter 

0.493 

0.077 

0.202 

0.031 

0.0236 

0.0013 

0.0179 

0.0005 

7 


Spring 

0.785 

0.373 

0.296 

0.135 

0.0236 

0.0021 

0.0191 

0.0009 

40 


Summer 

0.786 

0.312 

0.301 

0.129 

0.0242 

0.0028 

0.0191 

0.0015 

60 


Fall 

0.908 

0.453 

0.342 

0.159 

0.0234 

0.0016 

0.0187 

0.0007 

44 

25-60 m 

Winter 

0.267 

0.041 

0.121 

0.016 

0.0265 

0.0001 

0.0173 

0.0004 

3 


Spring 

0.331 

0.060 

0.127 

0.025 

0.0273 

0.0015 

0.0200 

0.0011 

20 


Summer 

0.376 

0.132 

0.135 

0.055 

0.0286 

0.0024 

0.0213 

0.0016 

55 


Fall 

0.319 

0.054 

0.127 

0.024 

0.0276 

0.0012 

0.0192 

0.0012 

19 

60-100 m 

Winter 

0.156 

0.061 

0.071 

0.027 

0.0306 

0.0021 

0.0178 

0.0010 

3 


Spring 

0.263 

0.038 

0.102 

0.016 

0.0284 

0.0014 

0.0199 

0.0010 

7 


Summer 

0.279 

0.142 

0.106 

0.055 

0.0311 

0.0026 

0.0209 

0.0019 

6 


Fall 

0.227 

0.023 

0.089 

0.015 

0.0289 

0.0006 

0.0195 

0.0009 

3 

100-1000 m 

Winter 

0.160 

0.037 

0.072 

0.018 

0.0298 

0.0023 

0.0185 

0.001 1 

5 


Spring 

0.178 

0.048 

0.071 

0.015 

0.0306 

0.0027 

0.0194 

0.001 1 

2 


Summer 

0.195 

0.124 

0.070 

0.051 

0.0338 

0.0036 

0.0223 

0.0030 

8 


Fall 

0.230 

0.110 

0.084 

0.039 

0.0305 

0.0025 

0.0204 

0.0006 

3 
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and Chi within the blue spectrum, which includes the Chi absorption 
peak near 443 nm, reduces the accuracy of satellite retrievals of Chi, 
especially for the heritage band-ratio algorithms (Carder et al., 1991 ; 
Siegel et al„ 2013; Siegel, Maritorena, Nelson, & Behrenfeld, 2005; 
Siegel, Maritorena, Nelson, Behrenfeld, & McClain, 2005). indeed, con- 
current changes in satellite CDOM and Chi using the standard Chi prod- 
uct (O'Reilly et al., 2000) provides false global-scale trends in Chi 
(Brown, Huot, Werdell, Gentili, & Claustre, 2008; Siegel et al., 2013). 
The capability to accurately retrieve (3 cdom(M provides the potential 
for improving satellite retrievals of Chi. 

Recent findings demonstrating the utility of the CDOM spectral slope 
at various wavelength ranges (Fichot & Benner, 2012; Helms et al., 
2008; Loiselle et al., 2009) have prompted greater interest in satellite re- 
trieval of S. For example, the ratio of S 2 75:295 and S 350: 4 oo is a very good 
indicator of DOM photochemical versus microbial degradation pro- 
cesses (Helms et al„ 2008). Satellite retrieval of S 30 o:6oo and Ocdom 
at a single wavelength within that spectral range allows computa- 
tion of Ocdom at any wavelength within the 300:600 nm spectral 
ranges using the CDOM spectral slope equation (Bricaud et al., 
1981). Furthermore, strong correlations between S275.295 and lignin 
phenol yields (Fichot & Benner, 2012) can be used to derive lignin 
phenol yields from satellite retrievals of S275.295 (Fichot et al., 2013) 
for subsequent estimation of vascular plant contributions to DOM in 
surface ocean waters. 

Ocean color satellite retrievals of a CDOM and S have many potential 
applications. Thus, the objectives of our work are (1) to develop bio- 
optical algorithms using field data from our study region within the 
northeast U.S. continental shelf, (2) validate these algorithms for 
SeaWiFS and MODIS-Aqua sensors, (3) examine what factors add 
significant uncertainty in satellite retrievals of o C dom and S, (4) de- 
scribe the distributions of o CD om and S from SeaWiFS and MOD1S re- 
trievals, and (5) evaluate whether radiometric measurements in the 
UV can improve the fidelity of o CD om and S algorithms. The results 
presented here are drawn from field sampling activities supported 
through multiple research grants in order to accumulate sufficient 
quantities of algorithm-independent coincident observations of satel- 
lite data and field observations. 

2. Methods 

2.J. Study area 

Measurements for this study were collected throughout the lower 
estuaries and continental margin of the northeastern U.S. between 
Cape Hatteras, North Carolina and southern Nova Scotia, Canada, 
which comprises the northeast U.S. continental shelf Large Marine 
Ecosystem (NES LME). For the purposes of this study and on the basis 


of ocean circulation, geography and bathymetry, the continental margin 
of the northeastern U.S. is subdivided into three regions (Fig. 1). The 
Middle Atlantic Bight (MAB) extends from Cape Hatteras, North 
Carolina to Cape Cod, Massachusetts. The Gulf of Maine (GoM) extends 
from Cape Cod to southern Nova Scotia, Canada. Georges Bank (GB) 
consists of the region to the southeast, east and northeast of Cape Cod 
(Fig. 1 ). The MAB is further sub-divided into the northern and southern 
MAB as shown in Fig. lc. 

2.2. Field sampling 

Multiple research cruises were conducted in the southern MAB from 
March 2005 to August 2007, lower Chesapeake Bay from May 2004 to 
August 2007, Gulf of Maine from April to June 2007, New York Bight 
(Hudson-Raritan Bay and adjacent coastal ocean) from May 2007 to 
May 2009, and across all continental shelf regions from August 2009 
to June 2011 (Fig. 1; Table 1). Seawater samples were collected with 
Go-Flo bottles at multiple depths per station. Additional samples were 
collected via peristaltic pumping from -0.5 m to 1 m depth to obtain 
coincident samples with in-water radiometry profiles as well as through 
the ships' clean seawater flow-through systems. 

2.3. Sample processing, storage, and analysis 

Seawater samples for measurement of CDOM absorbance spectra 
are collected, filtered, stored and analyzed as described in Mannino 
et al. (2008) and described here in abridged form. Samples are filtered 
through combusted Whatman GF/F glass fiber filters by gentle vacuum 
and collected directly into glass storage bottles. Samples are kept on ice 
for transport or in a refrigerator (4 to 8 °C) for longer-term storage. For 
analysis, CDOM samples are warmed to room temperature and filtered 
through either 0.2 pm Gelman Supor (polyethersulfone) filters or 
Whatman Nuclepore (polycarbonate) filters, which are pre-rinsed 
with ultrapure water (Milli-Q) and sample water. CDOM absorbance 
spectra are measured on a double-beam Caiy 100 Ultraviolet-Visible 
scanning spectrophotometer using Suprasil quartz cylindrical cells of 
100 mm pathlength with ultraviolet oxidized Milli-Q water as the 
blank and reference (Mannino et al., 2008; Mitchell, Kahru, Wieland, & 
Stramska, 2003). Scan settings for the spectrophotometer are as 
follows; 4 nm slit width, 250-800 nm wavelength range, 1 nm data in- 
terval, and 100 nm min -1 scan rate. The instrument noise for reference 
baseline spectral scans of air-to-air and Milli-Q water is evaluated daily 
to be within ± 0.0005 absorbance units. Typically, raw absorbance mea- 
surements are corrected with field filtration blanks of UV-oxidized 
Milli-Q and a null point value (Mitchell et al., 2000, 2003). Since the 
absorbance of filtration blanks and null point values are typically 
within the level of instrument noise, corrections are applied on a 


Table 3 

Coefficients and goodness of fit parameters for band-ratio bio-optical Qcdom algorithms developed from field observations of OcdomCR) and remote sensing reflectance (R rs (\)). R 2 is the 
non-linear regression coefficient. Sy.x represents the standard deviation of the residuals reported in units of the Y-axis variable, and n is the sample size. The minimum band ratio repre- 
sents the threshold R rs band ratio value at which the algorithm loses sensitivity below this value (at high a CDOM values). 


Exponential one-phase decay equation 

Qcdom275 

Qcdom3 55 

Qcdom3 80 

Gcdom412 

Gcdom443 

X = Ln[(Y- Bo) / B 2 ] / — B, 

B 0 = 0.2792 

0.2652 

0.2676 

0.2675 

0.2678 


Bj = 1.582 

5.534 

8.484 

13.74 

23.28 

Y = R rs (412) / R rs (547) 
X = Ocdom( X) 

B 2 = 21.95 

4.337 

4.054 

3.619 

3.406 

R 2 

0.880 

0.934 

0.944 

0.947 

0.928 

Sy.x 

0.099 

0.0838 

0.077 

0.0757 

0.082 

n 

140 

152 

152 

153 

150 

Minimum band ratio 

0.31 

0.295 

0.295 

0.295 

0.295 

X = Ln[(Y — B 0 ) / B 2 ] / — B, 

B 0 = 0.9686 

0.7723 

0.685 

0.7074 

0.7857 

Y = R„(412) / R rs (670) 

Bi = 2.302 

7.794 

9.522 

15.86 

31.79 

X = Ocdom( X) 

B 2 = 958.4 

92.44 

47.35 

43.85 

56.59 

R 2 

0.844 

0.943 

0.954 

0.946 

0.936 

Sy.x 

0.986 

0.817 

0.692 

0.605 

0.749 

n 

138 

146 

140 

137 

144 

Minimum band ratio 

1.29 

1.1 

1.1 

1.1 

1.1 
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limited set of samples. For the most recent sample collections 
(August 2009 to the present), a null point correction value is applied 
to all samples from the mean optical density within 695-700 nm. The 
naperian CDOM absorption coefficients are calculated from the fol- 
lowing equation: 

a c dqmW = 2.303A(\)/L (1) 

where A(\) is the absorbance of filtered seawater at a specific wave- 
length measured across pathlength L in meters. The CDOM spectral 
slope coefficients (S) are determined by fitting a single-exponential 
non-linear curve on multiple spectral ranges for each Ocdom wave- 
length range (275-295 and 300-600 nm; S 275.295 and S 300:B00 ) 
using the following expression: 

a(A) = a(\ 0 )e“ s(x " Xo) (2) 

where a(\) and a(\ 0 ) represent the absorption coefficients at wave- 
length A and reference wavelength X 0 . 

2.4. Apparent optical properties (AOP) 

Radiometry profiles were collected using either the Biospherical 
Profiler (BioPRO) or the Submersible Biospherical Optical Profiling 
System (SuBOPS). BioPRO is a compact rocket-shaped device, whereas 
the SuBOPS device uses a small kite-shaped backplane (Hooker, Lind, 
Morrow, & Brown, 2010). Both instruments were floated away to avoid 
deployment-platform perturbations and use the same 19-channel sen- 
sors spanning 320-875 nm with approximately 10 decades of dynamic 
range and a 12 Hz sampling rate. The adjustable buoyancy permits stable 
descents with vertical tilts to within 5°. The SuBOPS uses a hydrobaric 
buoyancy system and adjustable sensor orientation to achieve stable 
descents as slow as 5 cm s _1 with vertical tilts to within 2.5°. SuBOPS 
can achieve 1 cm (or less) vertical resolution in near-surface waters, 
because the compressible bladders in the hydrobaric buoyancy sys- 
tem allow the profiler to loiter close to the surface before reaching 
terminal velocity (usually set to 3-5 m). Both instruments possess 
19 micro-radiometers each with 10 nm-wide bands centered at the 
following wavelengths: 320, 340, 380, 395, 412, 443, 465, 490, 510, 
532, 555, 560, 625, 665, 670, 683, 710, 780, and 860 or 870 nm. 

The significance of acquiring useful data close to the sea surface is 
expressed directly in the processing scheme used to derive data prod- 
ucts from the light measurements. The processor used here is based 
on a well-established methodology (Smith & Baker, 1984) that was 
evaluated in an international round robin (Hooker et al., 2001) and 
shown to be capable of agreement at the 1% level when the processing 
options were as similar as possible. Complete details for the terms and 
dependencies are available in the Ocean Optics Protocols (hereafter, 
the Protocols), which initially adhered to the Joint Global Ocean Flux 
Study (JGOFS) sampling procedures (Joint Global Ocean Flux Study, 
1991 ) and defined the standards for NASA calibration and validation ac- 
tivities (Mueller & Austin, 1992). Over time, the Protocols were initially 
revised (Mueller & Austin, 1995), and then updated annually (Mueller, 
2000, 2002, 2003). 

The Protocols are detailed, so only a brief overview for obtaining 
data products from vertical profiles of upwelling radiance ( L u ) plus 
downward irradiance (£ d , respectively) are presented here. In-water ra- 
diometric quantities in physical units, 0 (i.e„ L u or £ d ), are normalized 
with respect to simultaneous measurements of the global solar irradi- 
ance, £ d (0 +, A, t), with t explicitly expressing the time dependence, 
according to 

j3(z, A, t 0 ) = 0(z, A, t)£„(0+, A, t 0 )/£ d (0+, z, t) (3) 

where /3(z, A, t 0 ) identifies the radiometric parameters as they would 
have been recorded at all depths z at the same time t 0 , and t 0 is generally 
chosen to coincide with the start of data acquisition. For simplicity, the 


variable t is omitted in the following text. In addition, any data collected 
when the vertical tilt of the profiler exceeds 5° are excluded from the 
ensuing analysis. 

After normalization and tilt filtering, a near-surface portion of 
£ d (z, A) centered at z 0 and having homogeneous optical properties 
(verified with temperature and attenuation parameters) extending 
from Zi = z 0 + A z and z 2 = z 0 — A z is established separately for 
the blue-green and red wavelengths; the ultraviolet (UV) is included 



a cooM(412) (nr 1 ) 



a CDOw(412) (nr 1 ) 


Fig. 2. Representation of band ratio algorithms from field measurements of the absorption 
coefficient of colored dissolved organic matter [a C DOM(412)] and remote sensing reflec- 
tance (Rrs) band ratios (a) R rs (490)/R re (555), (b) R rs (412)/R rs (555), and (c) R rs (412)/ 
R rs (665) for three coastal regions along the northeastern U.S., which include the Gulf of 
Maine, New York Bight (region between northern New Jersey and Long Island, New 
York) and Southern Middle Atlantic Bight (SMAB). The symbols represent the field obser- 
vations and the lines represent the predicted values produced from the non-linear model 
derived for each region or all regions combined. 
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in the interval most similar to the UV attenuation scales. Both intervals 
begin at the same shallowest depth, but the blue-green interval is 
allowed to extend deeper if the linearity in ln[L u (z, A)], as determined 
statistically, is thereby improved. The negative value of the slope of 
the regression yields the diffuse attenuation coefficient, K d (A), which 
is used to extrapolate the fitted portion of the E d profile through the 
near-surface layer to null depth, z = 0 _ . 

Fluctuations caused by surface waves and so-called lens effects pre- 
vent accurate measurements of £ d ( A) close to the surface. A value just 
below the surface (at null depth z = CP) can be compared to that mea- 
sured contemporaneously above the surface (atz = 0 + ) with a separate 
solar reference using 

£ d ((T,A) = 0.97£ d (o + ,A), (4) 

where the constant 0.97 represents the applicable air-sea transmit- 
tance, Fresnel reflectances, and the irradiance reflectance, and is deter- 
mined to an accuracy better than 1% for solar elevations above 30° 
and low-to-moderate wind speeds (Antoine, Hooker, Belanger, 
Matsuoka, & Babin, 2013). The distribution of E d measurements at any 
depth z influenced by wave focusing effects does not follow a Gaussian 
distribution, so linear fitting of E d in a near-surface layer is poorly 
constrained, especially if the number of samples is small. The applica- 
tion of Eq. (4) to the fitting process establishes a boundary condition 
or constraint for the fit. 

The appropriateness of the £ d extrapolation interval, initially 
established by z^ and z 2 , is evaluated by determining if Eq. (4) is 
satisfied to within approximately the uncertainty of the calibra- 
tions (a few percent); if not, z ^ and z 2 are re-determined while 
keeping the selected depths within the shallowest homogeneous 
layer possible until the disagreement is minimized (usually to 
within 5%). The linear decay of ln[/3(z, A)], for all light parameters 
in the chosen near-surface layer are then evaluated, and if linearity 
is acceptable, the entire process is repeated on a cast-by-cast basis. 
Subsurface primary quantities at null depth, /3(0 , A), are obtained 
from the slope and intercept given by the least-squares linear re- 
gression of ln[/3(z, A)] versus z within the extrapolation interval 
specified by z-i and z 2 . 

The water-leaving radiance is obtained directly from 



0.001 0.01 0.1 1 10 


[a,(412)-a,(443)] (nr') 



Fig. 3. Field measurement-based relationships of the total absorption coefficient [a f ] differ- 
ence of a t (412) and a f (443) with (a) a CDO M( 275), o C dom( 355), cicdom(412) and £Z CD om(443) 
and (b) CDOM spectral slope for two wavelength regions S 275295 and S 30 o: 6 oo- Qt{ 412) and 
a f (443) were computed from in situ R rs data using version 5 of the Quasi-Analytical 
Algorithm (QAA) from Lee et al., 2010, 2002. The symbols represent the field observations 
and the lines represent the predicted values produced from the non-linear models. 


L w (\) = 0.54L„(0 , A) (5) 

where the constant 0.54 accurately accounts for the partial reflection 
and transmission of the upwelled radiance through the sea surface, as 
confirmed by Mobley (1999). To account for the aforementioned de- 
pendence of L w on the solar flux, which is a function of atmospheric 
conditions and time of day, L w is normalized by the (average) global 


solar irradiance measured during the time interval corresponding to z x 
and z 2 : 

R rs (A)=L lv (A)/£ d (0 + ,A), (6) 

where R rs is the remote sensing reflectance. 


Table 4 

Coefficients and goodness of fit parameters for bio-optical CDOM algorithms developed from field observations of cicdom(^) and Quasi Analytical Algorithm (QAA) total absorption coef- 
ficient [a t (\)] derived from remote sensing reflectance. 


Equation 

Ocdom275 

Ocdom355 

Ocdom380 

Ocdom41 2 

Qcdom443 

S 275:295 

S 300:600 

QAA product-based 
Y = A * X B 

A = 10.08 

2.671 

1.771 

1.048 

0.617 

0.0158 

0.0162 


B = 0.4207 

0.575 

0.5829 

0.5797 

0.587 

-0.0737 

-0.1483 

X = a t (412) - a t (443) 

Y = Ocdom(^) or Scdom( ^ 1-^2) 








R 2 

0.884 

0.927 

0.929 

0.934 

0.928 

0.746 

0.881 

Sy.x 

0.664 

0.1467 

0.0962 

0.0547 

0.0339 

0.00088 

0.00144 

n 

150 

148 

148 

148 

148 

150 

150 

►< 

II 

> 

* 

A= 

1.236 

0.813 

0.499 

0.286 




B = 

0.600 

0.616 

0.646 

0.648 



X = a t (412) 

Y — °CDOm ( M 








R 2 


0.805 

0.811 

0.838 

0.826 



Sy.x 


0.198 

0.130 

0.0754 

0.0458 



n 


144 

144 

145 

145 
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2.5. CDOM algorithm development 

Spectral remote sensing reflectance derived from in situ radiometry 
profiles and o C domW were applied to develop multi-regional satellite 
algorithms to retrieve a CDOM and CDOM spectral slopes for the NES 
LME region. The station locations for measurements applied to derive 
algorithms of CDOM properties are shown in Fig. la (see Table 1 for 
additional details). A least squares approach was applied to field 
measurements to develop and evaluate multiple non-linear regression 
curve-fitting solutions as well as multiple linear regression analysis to 
retrieve a CD0M (^) from AOP-derived products ( R rs , I( d , inherent optical 
properties (IOP)). We regressed cicdom(X) and S with R rs band ratios, 
K d , and IOP products to deduce algorithms that are optimal for our 
study region. A series of absorption and backscatter products were de- 
rived from R rs data using version 5 of the Quasi-Analytical Algorithm 
(QAA v5; Lee et al., 2010, 2002). The QAA products applied for algorithm 
development include the total absorption coefficient (a f ) at 412 nm and 
the difference in a f (412) and a t (443), because the total absorption 
and backscatter products are considered to be the most robust products 
as opposed to the constituent products of phytoplankton (a p /,) and 
colored detritus absorption coefficients ( a CDM ). The QAA-derived 
[a f (412) — a f (443)] was applied to reduce the contribution of phyto- 
plankton absorption, because phytoplankton absorption is much higher 
at 443 nm than 412 nm, and this formulation would thus accentuate the 
CDOM absorption. The evaluation of the algorithm relationship is based 
on least squares statistics including the regression coefficient (R 2 ), the 
standard deviation of the data points from the non-linear regression 
curve fit (Sy.x) or the mean square error (MSE) and visual inspection 
of the in situ data and the modeled curve fit. Non-linear regression 
curve fitting and single linear regression analyses were conducted 
using Prism 5.0b (GraphPad Software, Inc.). The statistical software 
package SYSTAT 13 (Systat Software) was applied for the multiple linear 
regression analysis. 

2.6. Satellite validation of CDOM algorithms 

SeaWiFS merged local area coverage (MLAC) and MODIS-Aqua local 
area coverage (LAC) Level 1 files of the Northeast Atlantic are obtained 
from the NASA Ocean Color Web (http://oceancolor.gsfc.nasa.gov/). 
All files are processed to Level 2 using the standard defaults of SeaDAS 
6.4, and masks are applied to pixels with any of the followings flags: 
land, cloud or ice, stray light, sun glint, high top-of-atmosphere radiance, 
low normalized water-leaving radiance at 555 or 547 nm, or atmospheric 
correction failure, according to Bailey and Werdell (2006). Satellite 
validation data collected within 3 h of the in situ sampling are extracted 
in 3 x 3 native resolution pixel arrays centered on the field station loca- 
tions (Fig. lb). With the exception of the pixel array size extracted, all 
other satellite validation protocols are in accordance with those 
described in Bailey and Werdell (2006). For all the CDOM validation 
matchups, the average distance between the field station location and 
the center point of the center pixel in the 3x3 validation pixel array 
is 0.66 km. The validation of algorithms is conducted with a completely 
independent dataset from the algorithm development dataset, i.e. data 
used in algorithm development are excluded from the validation analy- 
sis (Table 1; Fig. t). 

The evaluation of algorithm performance was based on statistical pa- 
rameters comparing the satellite-derived retrievals of products with the 
field measurements, which are referred to here as validation matchups. 
The filtered means from the satellite 3x3 pixel arrays were computed 
for comparison with field measurements: 

Filtered Mean = Ej(1.5 * cr— X)<Xj<(1.5*cr + X)]/N (7) 

where X is the unfiltered mean value, a is the standard deviation of the 
unfiltered values, x, are the individual values for each 3x3 pixel 
array, and N is the number of values from each 3x3 pixel array that 


falls within ± 1.5 *cj (Bailey & Werdell, 2006). The coefficient of varia- 
tion (CV) is computed from the filtered mean and standard deviation 
of values that pass the filtered criteria. Validation sites that exceed the 
mean filtered CV by >0.15 are excluded from the validation analysis. A 
minimum of 5 pixels from the 3x3 pixel arrays is required for inclusion 
of the matchup in the validation analysis. 
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Fig. 4. Comparison of field measurements of (a) Ocdom( 412), (b) S275.295 and (c) S 30 o:6oo 
with predicted values from multiple linear regression (MLR) algorithm of R„(443) and 
R rs (547) for the red circles labeled MLR or with ultraviolet and other bands (MLR-UV) 
shown as blue triangles, which include R rs (380), R„(412), R rs (443), R„(490), Rrs(532), 
R„(547) and R rs (665) for plots (a) and (c) and R„(380), R rs (412), R rs (443), R rs (490) and 
R rs (547) for plot (b). (For interpretation of the references to color in this figure legend, 
the reader is referred to the web version of this article.) 
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Table 5 

Coefficients and goodness of fit parameter values from multiple linear regression-based (MLR) bio-optical a CDOM algorithms developed from field observations of £1 cdom(A) and remote 
sensing reflectance (Rj^K)). Y = Ho + B i X] + B2X2 -B B3X3 -B B4X4 -B B5X5 -B B^Xii -B B7X7; Xj 7 = Ln[R rs ( X 7 . . . X 7 ) ] ; Y = Ln[ticDOM(X)] or Ln [Sqxjm( X 1 ! X2 ) ]■ Adj. R 2 is the adjusted 
R 2 , and MSE is the mean square error of the multiple linear regression analysis. 


Multiple linear regression 

Qcdom275 

Gcdom355 

Gcdom380 

Qcdom4 12 

Qcdom443 

S 275:295 

5300:600 

MODIS bands 

B 0 = constant 

0.464 

-1.960 

-2.507 

-3.070 

-3.664 

-3.258 

-3.640 

B,: R rs 443 

-0.769 

-1.208 

-1.261 

-1.285 

-1.291 

0.336 

0.186 

B 2 : R rs 547 

0.692 

1.049 

1.088 

1.107 

1.105 

-0.279 

-0.146 

Adj. R 2 

0.680 

0.764 

0.778 

0.790 

0.775 

0.692 

0.748 

MSE 

0.057 

0.090 

0.090 

0.084 

0.092 

0.0096 

0.0021 

n 

155 

155 

155 

154 

154 

155 

155 

SeaWiFS bands 

B 0 = constant 

0.643 

-1.692 

-2.227 

-2.784 

-3.379 

-3.325 

-3.679 

Bi: R rs 443 

-0.682 

-1.076 

-1.124 

-1.146 

-1.1513 

0.300 

0.168 

B 2 : R rs 555 

0.630 

0.954 

0.990 

1.008 

1.006 

-0.252 

-0.134 

Adj. R 2 

0.670 

0.751 

0.765 

0.777 

0.762 

0.674 

0.740 

MSE 

0.059 

0.095 

0.095 

0.089 

0.097 

0.010 

0.0022 

n 

155 

155 

155 

154 

154 

155 

155 

uv and MODIS bands 

B 0 = constant 

0.4467 

-2.092 

-2.677 

-3.176 

-3.819 

-3.085 

-3.601 

Bi ! R rs 380 

-0.5358 

-0.604 

-0.598 

-0.530 

-0.557 

0.0747 


B 2 : R rs 412 

-0.4819 

-1.265 

-1.319 

-1.423 

-1.510 

0.824 

0.321 

B 3 : R rs 443 

1.4978 

2.575 

2.667 

2.714 

3.000 

-1.416 

-0.449 

B 4 ; R rs 490 

-1.771 

-2.479 

-2.502 

-2.513 

-2.776 

1.075 

0.335 

B 5 : R rs 532 

1.475 

1.309 

0.803 

0.681 

0.715 


0.222 

B 6 : R rs 547 
B 7 : R rs 665 

-0.4864 

0.1726 

0.215 

0.491 

0.190 

0.602 

0.200 

0.659 

0.185 

-0.459 

-0.375 

Adj. R 2 

0.834 

0.916 

0.926 

0.931 

0.924 

0.910 

0.847 

MSE 

0.030 

0.032 

0.030 

0.027 

0.031 

0.0028 

0.0013 

n 

155 

155 

155 

154 

154 

155 

155 


The statistical parameters applied for algorithm validation include 
the mean and standard deviation of the absolute percent difference 
(APD), root mean square error (RMSE), median ratio of computed fil- 
tered mean satellite value (C a i g ) to field measurement (C in situ), the 
semi-interquartile range (SIQR) and the R 2 and slope values from linear 
regression analyses of the validation matchups for each satellite sensor 
(Bailey & Werdell, 2006). 


Mean APD(X) = [S | (c alg - C in situ ) |/C in situ ] * 100/N (8) 


RMSE = 


1/N*l£(c alg -C insitu ) 


2 ] 1 / 2 


( 9 ) 


Table 6 

Coefficients and goodness of fit values for bio-optical Ocdom algorithms developed from 
field observations of a CD 0 M(M and light attenuation coefficients (K d (\)). 


Equation 

Qcdom 355 

Qcdo/v?380 

Ocdom 412 

Qcdom443 

-< 

II 

> 

>< 

A = 0.5097 

0.3307 

0.1979 

0.1145 

X = 7C d (340) 

B = 0.9321 

0.9431 

0.936 

0.9449 

Y — Ocdom(^) 

R 2 

0.952 

0.953 

0.955 

0.947 

Sy.x 

0.1153 

0.07594 

0.0437 

0.0283 

n 

148 

148 

148 

149 

-< 

II 

> 

>< 

A = 0.8325 

0.5409 

0.3207 

0.187 

X = K,i(380) 

B = 0.7928 

0.8001 

0.7961 

0.8017 

y — Qcdom(^) 

R 2 

0.884 

0.8911 

0.905 

0.892 

Sy.x 

0.1818 

0.1157 

0.0635 

0.0402 

n 

152 

151 

150 

151 

Y = A . X B 

A = 1.021 

0.6680 

0.4006 

0.2311 

X = K„(412) 

B = 0.7076 

0.7165 

0.7141 

0.72 

y — Qcdom(^) 

R 2 

0.818 

0.8217 

0.8472 

0.8244 

Sy.x 

0.2151 

0.1416 

0.0770 

0.0490 

n 

150 

150 

149 

150 




Fig. 5. Relationships of light attenuation coefficient at 340 nm (a) /C d (340) and at 412 nm 
(b) K d { 412) with Ocdom( 3 55) and a CD 0 M(4 12) from field observations. 
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* Bias = [(l/N * E (C alg -C in Situ ))/Mean (Q n situ )] * 100 (10) 

SIQR=(Q3-Qi)/2 (11) 

where Q 3 and Ch represent the 75th and 25th percentiles for the ratios 
of the satellite-derived values to the field measurement. 

3. Results 

3.1. Field distributions of CDOM properties 

The spatial and seasonal distributions of a CDOM , ^ 275.295 and S 300: soo 
surface measurements collected from our cruises throughout the 
NES LME are presented in Table 2 (see Fig. lc for a map of the sub- 
regions). These results show expected trends of decreasing o C dom and 
increasing S 275.295 and S 300:600 between the shore and continental 
slope for the Gulf of Maine and both MAB regions. The shallower areas 
of Georges Bank, a strong tidally mixed region, have higher a CD0M and 
lower S 27 5:295 and S 30 o: 6 oo than the deeper portions of Georges Bank. In 
the northern MAB (NMAB), a CDOM on the inner- and mid-shelf peaked 
in spring while 3275.295 was lowest, which coincides with peak river dis- 
charge. However, in the southern MAB (SMAB), a CDOM was higher on av- 
erage on the inner shelf in the fall and on the mid-shelf during summer. 
Note the high standard deviations for these values reaching -50% 


indicating high spatial and temporal variability in Ocdom- Because 
some seasons and sub-regions were either not sampled or represented 
by too few stations and collected over 7 years, it is not possible to draw 
too many conclusions from these results especially for the outer shelf 
and continental slope areas and Gulf of Maine (Table 2), hence the 
utility of satellite observations. 

32. Algorithm development 

3.2.1. Band-ratio OcdomM algorithms 

Multiple R rs band-ratios yield strong relationships with a CDOM {\ ) for 
our northeastern U.S. coastal study region (Table 3). Prior work on 
a CD oM algorithm development and validation for the southern MAB 
recommended a band-ratio algorithm of R rs (490)/R rs (555) and an 
exponential one-phase decay non-linear model (Mannino et al., 2008). 
However, as new field observations from the Gulf of Maine, Hudson 
estuary plume and New York Bight were added to the southern MAB 
dataset, it became clear that this band-ratio algorithm would not be 
appropriate across multiple regions (Fig. 2a). Applying the same 
non-linear model with a different band ratio, R rs (412)/R rs (555), 
yielded a very good fit to the field observations for q cdom ( 412) 
(Fig. 2b; R 2 = 0.95; Sy.x = 0.0806) and other o C domM (Table Al). 
Other band ratios such as R rs (412)/R rs (547), R rs (412)/R rs (670), and 
R rs (412)/R rs (665) (Fig. 2c) also yielded strong relationships for a broad 
range of q C dom(^) (Table 2; Table A2). Although these algorithms have 


Table 7 

Satellite validation statistics of CDOM algorithms for SeaWiFS and MODIS-Aqua. Values and algorithms shown in bold represent algorithms with the best validation results, sd is the stan- 
dard deviation of the MAPD. SIQR is the semi-interquartile range. MLR_Fichot refers to an MLR algorithm published by Fichot et al. (2013). 


Product 

Algorithm 

MAPD 

±sd 

RMSE 

Median ratio 

SIQR 

%Bias 

Slope 

R 2 

n 

SeaWiFS 

Qcdom 275 

Exp_R„(412)/R rs (670) 

26.2 

24.9 

1.76 

0.94 

0.182 

-24.6 

0.098 

0.19 

33 


QAA-based 

43.4 

26.1 

2.03 

0.91 

0.417 

-19.7 

0.127 

0.04 

43 


MLR 

27.3 

20.6 

1.49 

0.95 

0.225 

- 18.4 

0.395 

0.53 

42 

Qcdom355 

Exp_R rs (412)/R rs (670) 

28.1 

29.4 

0.329 

0.96 

0.189 

-25.9 

0.133 

0.21 

33 


QAA-based 

51.5 

33.2 

0.410 

0.86 

0.464 

- 21.2 

0.168 

0.06 

43 


MLR 

29.0 

21.9 

0.305 

0.78 

0.215 

-29.0 

0.541 

0.59 

42 

Qcdom 3 80 

Exp_R„(412)/R rs (670) 

27.5 

27.4 

0.200 

0.94 

0.218 

-26.4 

0.166 

0.21 

33 


QAA-based 

51.4 

32.7 

0.262 

0.80 

0.450 

- 20.8 

0.179 

0.06 

43 


MLR 

29.1 

21.0 

0.196 

0.75 

0.203 

-30.5 

0.567 

0.61 

42 

q cdom4 12 

Exp_R rs (412)/R rs (670) 

27.7 

25.7 

0.114 

0.90 

0.192 

-27.1 

0.179 

0.22 

33 


QAA-based 

50.6 

32.9 

0.152 

0.84 

0.457 

- 20.1 

0.188 

0.064 

43 


MLR 

29.6 

19.0 

0.113 

0.75 

0.182 

-31.0 

0.593 

0.62 

42 

a CDOM ^ 43 

Exp_R rs (412)/R rs (670) 

28.4 

25.6 

0.0687 

0.86 

0.163 

-29.5 

0.165 

0.23 

33 


QAA-based 

51.0 

33.5 

0.0877 

0.86 

0.472 

- 20.1 

0.192 

0.07 

43 


MLR 

30.0 

19.0 

0.0650 

0.76 

0.192 

-31.0 

0.58 

0.62 

42 

$275:295 

QAA-based 

18.9 

36.1 

0.0079 

1.03 

0.097 

9.1 

-0.174 

0.01 

43 


MLR 

9.9 

10.4 

0.00318 

1.06 

0.042 

7.5 

0.545 

0.44 

42 

$ 300:600 

QAA-based 

8.3 

9.1 

0.00240 

1.03 

0.053 

6.0 

0.642 

0.19 

43 


MLR 

9.1 

6.0 

0.00211 

1.07 

0.045 

7.2 

0.264 

0.10 

42 

MODIS 

Qcdom 275 

Exp_R„(412)/R rs (547) 

23.4 

17.1 

0.952 

1.09 

0.182 

- 1.9 

0.270 

0.37 

53 


QAA-based 

53.1 

41.3 

1.504 

1.37 

0.393 

25.6 

0.738 

0.47 

55 


MLR 

24.3 

18.1 

1.086 

1.01 

0.214 

-8.2 

0.532 

0.58 

61 

Qcdom355 

Exp_R rs (412)/R rs (547) 

27.1 

20.6 

0.186 

1.01 

0.251 

- 4.0 

0.360 

0.39 

53 


QAA-based 

69.0 

56.4 

0.374 

1.48 

0.469 

36.5 

1.05 

0.52 

55 


MLR 

23.9 

19.1 

0.212 

0.87 

0.150 

-17.4 

0.626 

0.76 

61 

Qcdom 3 80 

Exp_R„(412)/R rs (547) 

25.9 

19.5 

0.118 

0.99 

0.225 

-5.4 

0.379 

0.40 

53 


QAA-based 

69.2 

56.3 

0.250 

1.52 

0.466 

37.9 

1.10 

0.52 

55 


MLR 

23.6 

18.4 

0.134 

0.85 

0.163 

-18.5 

0.658 

0.76 

61 

a CDOM^ 12 

Exp_R„(412)/R rs (547) 

25.1 

18.7 

0.0674 

0.95 

0.194 

-7.5 

0.415 

0.42 

53 


QAA-based 

66.7 

55.6 

0.149 

1.53 

0.438 

37.7 

1.13 

0.51 

55 


MLR 

22.7 

17.0 

0.0760 

0.85 

0.147 

-19.1 

0.702 

0.77 

61 

q cdo/w443 

Exp_R rs (412)/R rs (547) 

24.7 

18.1 

0.0382 

0.95 

0.186 

-7.5 

0.434 

0.43 

53 


QAA-based 

67.4 

54.8 

0.0881 

1.52 

0.423 

39.3 

1.17 

0.52 

55 


MLR 

20.4 

14.5 

0.0430 

0.86 

0.138 

-18.6 

0.706 

0.77 

49 

$ 275:295 

QAA-based 

11.4 

10.3 

0.00390 

0.92 

0.072 

-4.6 

0.723 

0.19 

55 


MLR 

8.7 

6.4 

0.00264 

1.06 

0.044 

6.7 

1.02 

0.64 

59 


MLR_Fichot 

7.2 

5.8 

0.00225 

1.02 

0.054 

2.0 

1.05 

0.54 

52 

$ 300:600 

QAA-based 

6.3 

7.2 

0.00180 

1.01 

0.035 

1.4 

0.261 

0.04 

55 


MLR 

9.7 

4.2 

0.00197 

1.10 

0.036 

8.2 

0.649 

0.32 

59 
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very high sensitivity for the lower a CDOM values measured within our 
study region, they have little sensitivity at higher o C dom values typical 
of nearshore waters including the western Gulf of Maine, Hudson River 
Estuary, Chesapeake Bay and Delaware Bay. For instance, there is only a 
small change in the R rs (412)/R rs (555) value for o CDO m(412) values 
between -0.3 and 1 m ' 1 of the non-linear model. A cubic polynomial 
fit of the natural log of a CD0M (\) with R rs (412)/R re (555) yields similar re- 
sults as the exponential one-phase decay model (results not shown). 
Band ratio algorithms where the blue band is replaced with an ultraviolet 
band (e.g„ 340, 380 or 395 nm) yield similarly good relationships (data 
not shown), but also poor sensitivity at high a CDOM values as the blue/ 
green and blue/red band ratios. 

3.2.2. IOP-based CDOM algorithms 

The most promising set of algorithms relate a CDOM W to absorption 
products derived from QAA v5 (Lee et al„ 2002, 2010). There is a strong 
positive relationship between a CD0M ( 412) and the total absorption 
coefficient at 412 nm [a t (412)] derived using QAA from in situ R rs , 
but a great deal of scatter at higher a CDO jvi(412) values (Table 4; 
R 2 = 0.84; Sy.x = 0.0754). In contrast, the non-linear model fit of 
[a f (412) — a f (443)] versus a CD0M (\ ) yielded a strong relationship 
across nearly three orders of magnitude of o C dom (Table 4; Fig. 3a). 
Although results for only five a CDO M wavelengths are shown, the ap- 
proach is applicable across a broad range of wavelengths such as from 
250 to 450 nm. This approach also yields strong relationships between 
[a f (412) — a t (443)] and CDOM spectral slopes, S 27 s : 2 g5 and S 30 o: 6 oo 
(Fig. 3b; Table 4; R 2 = 0.75; Sy.x = 0.00088 and R 2 = 0.88; 
Sy.x = 0.0014, respectively). The S 2 75:295 relationship shows substan- 
tial scatter at high S 2 7 5:2g5 values and thus may not be appropriate for 
[a f (412) — a t (443)] greater than 0.1 m -1 . The S 30 o : boo relationship 
loses sensitivity below 0.02 nm -1 . These algorithms are referred to as 
QAA-based algorithms throughout the remainder of this document. 

3.2.3. Multiple linear regression CDOM algorithms 

Bio-optical algorithms for several CDOM products were developed 
by applying multiple linear regression (MLR) analysis from field mea- 
surements of multiple R rs bands to retrieve o C dom(A). S 275.295 and 
S 30 o:6oo- Initially, a stepwise backward statistical approach was applied 
using field-derived R rs bands close to the band centers found on the 
SeaWiFS and MODIS-Aqua satellite sensors (443, 490, 510 or 531, 547 
or 555, and 665 or 670 nm). The 412 nm band was not applied due to 
concerns of the high uncertainties in satellite retrievals of this particular 
band (e.g., Goyens, Jamet, & Schroeder, 2013). The MLR equation takes 
the following form: 


show a significant improvement from the MLR analyses that did not 
include the UV and 412 nm band (e.g., adjusted R 2 = 0.93 and 
MSE = 0.027 fora CDOM (412); adjusted R 2 = 0.91 and MSE = 0.0028 
for S 27 5 . 295 ; adjusted R 2 = 0.85 and MSE = 0.0013 for S 30 o: 6 oo', Fig. 4; 
Table 5). Plots of the field measurements versus the UV-MLR retrieved 
values demonstrate improvements in Qcdom(412) retrievals (and for 
S275.295 and S 300:600 ) across the full range of a CD0M ( 412) measurements 
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Y = B 0 + B,X, + B 2 X 2 + B 3 X 3 + ... B n X n (12) 

where X l n = Ln[R rs (X 1 ...\ n )],B 0 = constant, B^.Bn as coefficients for 
the R rs bands, and Y = Ln[a CDOM (\) ] or Ln(S). The validation statistics of 
the satellite-retrieved values and independent field data for the step- 
wise backward algorithms requiring three to five R rs bands were not 
as good as the guided MLR analysis using two R rs bands (443 and 547 
or 555 nm) (results not shown for the stepwise backward algorithms). 
The two-band MLR algorithm curve fitting results are encouraging (ad- 
justed R 2 = 0.79 and MSE = 0.084 for a CDOM { 412); adjusted R 2 = 0.75 
and MSE = 0.0022 for S 300:60 o', Fig. 4; Table 5), though not as strong as 
the regression results from the band-ratio and QAA-based algorithms. 

3.2.4. CDOM algorithms that employ ultraviolet radiometry 

Ultraviolet radiometric measurements were incorporated into 
CDOM bio-optical algorithms to determine whether UV information 
would improve the accuracy of CDOM algorithms. Thus, the stepwise 
backward MLR analysis was applied to include several more bands in 
the blue and UV spectral range including 380, 395 and 412 nm plus 
the 5 other bands (443, 490, 532, 547 and 665 nm) discussed previously 
(UV-MLR). The goodness of fit statistics for these UV-MLR analyses 



Fig. 6. Validation matchups of satellite-derived czcdom(4 12) products from SeaWiFS and 
MODIS-Aqua and field observations for (a) MLR, (b) QAA-based, and (c) exponential 
decay band ratios that applies the R „ band ratios R rs (412)/R re (547) for MODIS-Aqua and 
R re (412)/R re (670) for SeaWiFS (see Table 7 for validation statistics). 
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compared to the MLR retrievals, especially at the higher end of the 
range (Fig. 4). The MLR retrieved values shown in Fig. 4 were not ob- 
tained from an independent dataset from the data used to develop the 
MLR algorithms, but are included to demonstrate the potential of the 
MLR algorithms for both existing satellite sensors and planned future 
sensors. Since planned future sensors such as PACE and GEO-CAPE 
are not expected to provide R rs below -350 nm, the results for the 
UV-MLR with the 340 nm band are not shown here but are not signifi- 
cantly different than the results shown in Table 5. 

Algorithms for Ocdom(^) were also developed through non-linear re- 
gression with the light attenuation coefficients, I< d (\), for blue (41 2 nm) 
and UV bands (340 and 380 nm), which were computed from in-water 
radiometry profiles. The relationships with K d (340) yielded the 
strongest relationships with Ocdom(^) (for o C dom( 412) R 2 = 0.96 and 
Sy.x = 0.044), while the K d (412) relationship was weakest (R 2 = 0.85 
and Sy.x = 0.077; Table 6; Fig. 5). 


3.3. Satellite validation of a CDOM (\) and S algorithms 

Algorithm performance is evaluated with statistical parameters cal- 
culated from the filtered mean satellite data ( R rs and a f ) and field mea- 
surements (Table 7). For both SeaWiFS and MODIS-Aqua, we validated 
multiple algorithms each for o C dom( 275), q C dom(355), Qcdom( 380), 
Ocdom{ 412), and Qcdom( 443) and for other wavelengths (254, 295, 
300, 360 nm) for which data is not shown because results were similar. 
K d and UV-based MLR algorithms could not be validated due to the lack 
of necessary satellite products for Ocdom(A-) retrieval (K d or UV R rs 
bands). 



in situ S 275:295 (nm- 1 ) 



Fig. 7. Validation matchups of satellite-derived CDOM spectral slope products from SeaWiFS 
S 3 00 : 600 . and (d) QAA-based S 300:600 products (see Table 7 for validation statistics). 


The validation statistics indicate that the QAA-based algorithm did 
not perform as well for satellite retrieval of Qcdom(^) as compared to 
the MLR and two band-ratio algorithms. The mean absolute percent 
difference (MAPD) and RMSE were consistently greater for the 
QAA-based algorithms for both SeaWiFS and MODIS-Aqua, ranging 
from 43.4% ± 26.1 and 2.03 m _1 for Qcdom(275) to 51.4% ± 32.7 
and 0.26 m -1 for a CD 0 M( 380) for SeaWiFS, and from 53.1% ± 41.3 
and 1.50 m _1 for o C dom( 275) to 69.2% ± 56.3 and 0.25 m' 1 for 
Qcdom( 380) for MODIS (Table 7). For the QAA-based a CD0M (\ ) algo- 
rithms, both the MAPD and percent bias were greater for the MODIS 
validation than SeaWiFS. 

Both the MLR and band-ratio algorithms performed well in retrieval 
of o CDOM (\). Of the two band ratios, the Exp_R rs (412)/R rs (670) per- 
formed better for SeaWiFS retrieval of OcdomM, and Exp_R rs (412)/ 
R rs (547) performed better for MODIS (Table 7; Table A2). The preferable 
validation metrics of Exp_R rs (412)/R rs (670) for SeaWiFS are due to the 
exclusion of data values for which the SeaWiFS R rs ( 670) 3 x 3 pixel 
array values did not meet the filtering criteria. The MLR algorithms per- 
formed better than the others when all the validation metrics are taken 
into account, particularly the slope and regression values from the linear 
regression of field measurements versus satellite-retrieved values 
(Fig. 6; Table 7). On average for all five OcdomW MODIS retrievals, the 
MLR algorithm for MODIS yielded an MAPD of 23.0% ± 17.4, median 
ratio of 0.89, S1QR of 0.162, bias of — 16.4%, slope of 0.64 and R 2 of 
0.73. For SeaWiFS, the average values for the validation metrics were 
MAPD of 29.0% ± 20.3, median ratio of 0.80, SIQR of 0.203, bias 
of — 28%, slope of 0.54 and R 2 of 0.59. The validation results for MODIS 
OcdomW were generally better than for SeaWiFS, especially on bias, 
RMSE, slope and median ratio (Table 7; Fig. 6). 




in situ S m . 600 (nm 1 ) 

and MODIS-Aqua and field observations for (a) MLR S275295, (b) QAA-based S275295, (c) MLR 
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Table 8 

Validation of SeaWiFS and MODIS-Aqua remote sensing reflectance ( R rs ) products. 


Product 

In situ band center (nm) 

MAPD 

±sd 

RMSE 

Median ratio 

SIQR 

%Bias 

Slope 

R 2 

n 

SeaWiFS 

R„(412) 

412 

101.9 

86.7 

0.00254 

2.12 

0.723 

99.5 

0.992 

0.19 

18 

Rrs(443) 

443 

70.5 

53.1 

0.00203 

1.72 

0.520 

68.0 

1.02 

0.24 

19 

R rs (490) 

490 

36.5 

21.2 

0.00134 

1.21 

0.337 

15.5 

0.672 

0.22 

22 

R rs (510) 

510 

30.2 

18.1 

0.00110 

1.08 

0.279 

8.42 

0.618 

0.34 

25 

Rrs( 555) 

555 

29.7 

23.1 

0.00118 

1.03 

0.245 

-1.65 

0.583 

0.39 

28 

R„(670) 

670 

60.9 

71.3 

0.000343 

1.10 

0.435 

5.30 

0.732 

0.78 

23 

R rs (412)/R rs (555) 


57.7 

54.3 

0.437 

1.42 

0.222 

50.1 

1.02 

0.85 

18 

R rs (412)/R rs (670) 


44.4 

59.2 

2.35 

1.09 

0.227 

0.18 

0.700 

0.89 

18 

R rs (490)/R rs (555) 


11.9 

12.0 

0.136 

1.00 

0.105 

-1.25 

0.838 

0.94 

22 

MODIS-Aqua 

R„(412) 

412 

41.3 

26.5 

0.00113 

0.77 

0.233 

-3.0 

1.24 

0.47 

10 

R|-s (443) 

443 

26.6 

17.3 

0.000935 

0.85 

0.202 

-2.3 

1.56 

0.59 

11 

Rrs(469) 

465 

19.3 

10.1 

0.000799 

0.82 

0.150 

-2.7 

1.38 

0.69 

11 

Rrs(488) 

490 

19.9 

11.1 

0.000930 

0.77 

0.133 

-14.4 

0.951 

0.60 

12 

R rs (531) 

532 

14.5 

10.4 

0.000865 

0.86 

0.087 

-13.8 

0.817 

0.74 

13 

R„(547) 

Modeled 

15.3 

10.7 

0.000940 

0.85 

0.084 

-15.8 

0.764 

0.74 

13 

R„(555) 

555 

23.2 

12.8 

0.00114 

0.78 

0.0963 

-22.4 

0.761 

0.75 

14 

R rs (667) 

665 

31.2 

12.1 

0.000224 

0.67 

0.0993 

-21.0 

1.07 

0.81 

13 

R rs (667) 

670 

38.9 

11.2 

0.000293 

0.59 

0.0779 

-26.4 

1.04 

0.69 

13 

R rs (678) 

683 

38.4 

12.6 

0.000318 

0.65 

0.116 

-24.4 

0.891 

0.63 

13 

R rs (412)/R rs (547) 

As above 

36.6 

22.0 

0.316 

0.92 

0.293 

9.8 

1.30 

0.92 

9 

R„(412)/R rs (667) 

665 

42.9 

63.3 

5.667 

0.96 

0.266 

54.8 

1.55 

0.93 

9 

R„(488)/R rs (547) 

As above 

12.0 

7.6 

0.210 

0.93 

0.088 

1.2 

1.49 

0.94 

12 


Validation statistics for the CDOM spectral slope algorithms were 
quite good for both sensors and algorithms evaluated. Excluding the 
QAA-based S 275.295 algorithm, the validation statistics for the other 
S 275:295 and S 300:6 oo algorithms averaged 8.5% ± 7 for MAPD, 0.0023 


for RMSE, 1.05 for the median ratio, 0.044 for SIQR, 5.6% for bias, 0.63 
for slope and 0.32 for the R 2 (Table 7). For SeaWiFS, the QAA-based 
$ 300:600 algorithm performed slightly better statistically than the MLR 
algorithm, but the opposite was the case for MODIS (Fig. 7; Table 7). 


Table 9 

Spearman rank correlation values of SeaWiFS and MODIS-Aqua R rs validation residuals versus various SeaWiFS and MODIS-Aqua parameters (atmospheric, sensor geometry, pixel size, 
etc.). P-value adjusted for number of comparisons, p = 0.05/13 = 0.00385. 


SeaWiFS 

solz 

senz 

N0 2 tropo 

Epsilon 

AOT 490 

AOT 555 

AOT 865 

Angstrom 

Time Diff 

Pixel Size Y 

Pixel Size X 

Pixel Area 

R rs (670) 

R rs (412) 

-0.721* 

0.746* 

ns 

ns 

0.785* 

0.783* 

0.747* 

ns 

ns 

ns 

0.708* 

0.681* 

ns 

R rs (443) 

-0.711* 

0.670* 

ns 

ns 

0.806* 

0.802* 

0.767* 

ns 

ns 

ns 

0.638* 

0.61 1 A 

ns 

Rrs(490) 

-0.645* 

ns 

— 0.569 A 

ns 

0.705* 

0.701* 

0.668* 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (510) 

-0.586* 

ns 

-0.560* 

ns 

0.679* 

0.675* 

0.625* 

ns 

ns 

ns 

ns 

ns 

ns 

Rrs(555) 

-0.651* 

ns 

-0.623* 

ns 

0.761* 

0.757* 

0.714* 

ns 

ns 

ns 

ns 

ns 

ns 

Rrs(670) 

-0.596* 

ns 

-0.645* 

ns 

0.692* 

0.691* 

0.642* 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (412)/R rs (555) 

— 0.646 A 

0.637 A 

ns 

ns 

0.730* 

0.728* 

0.726* 

ns 

ns 

ns 

0.608 A 

0.584 A 

ns 

R rs (412)/R„(670) 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (490)/R rs (555) 

ns 

— 0.425 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

— 0.437 A 

0.537 A 

Mean 

34.9 

35.7 

4.89E + 15 

1.126 

0.220 

0.183 

0.091 

1.578 

1.4 

1.2 

2.1 

2.5 

0.00099 

sd 

13.9 

10.1 

3.43E + 15 

0.047 

0.169 

0.139 

0.065 

0.439 

0.78 

0.01 

0.69 

0.80 

0.00060 

n 

28 

28 

28 

28 

28 

28 

28 

27 

28 

28 

28 

28 

23 


MODIS-Aqua 

solz 

senz 

N0 2 tropo 

Epsilon 

AOT 488 

AOT 547 

AOT 869 

Angstrom 

Time Diff 

Pixel Size Y 

Pixel Size X 

Pixel Area 

R„(645) 

Rrs(412) 

ns 

ns 

— 0.783 A 

0.750 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (443) 

ns 

ns 

— 0.645 A 

0.691 A 

ns 

ns 

ns 

ns 

0.673 A 

ns 

ns 

ns 

ns 

Rrs(469) 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (488) 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (531) 

ns 

ns 

— 0.626 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

Rrs(547) 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (555) 

ns 

ns 

— 0.564 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (667) 

ns 

ns 

— 0.692 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (667) 

ns 

ns 

-0.850* 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (678) 

— 0.537 A 

ns 

-0.749* 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (412)/R rs (547) 

ns 

ns 

-0.883* 

0.833 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (412)/R„(667) 

ns 

ns 

— 0.717 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (412)/R„(555) 

ns 

ns 

-0.883* 

0.833 A 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

R rs (488)/R rs (547) 

ns 

ns 

-0.776* 

0.769* 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

ns 

Mean 

35.0 

45.4 

6.69 E + 15 

1.15 

0.16 

0.141 

0.079 

1.33 

1.18 

1.42 

2.70 

4.01 

7.47 E- 04 

sd 

12.7 

9.5 

3.18 E + 15 

0.10 

0.10 

0.085 

0.047 

0.55 

0.84 

0.21 

0.83 

1.87 

5.09 E- 04 

n 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

15 

13 


solz = solar zenith angle, sens = sensor zenith angle; N0 2 tropo = tropospheric column N0 2 concentration; epsilon = ratio of aerosol reflectance for two NIR bands applied in aerosol 
model selection; AOT = aerosol optical thickness at wavelengths (nm) indicated; Angstrom = angstrom exponent of the aerosol model applied in LI to L2 SeaDAS processing; Time 
Diff = time difference in hours between satellite observation and in situ sample collection; Pixel Size Y = pixel dimension from north to south in km; Pixel Size X = pixel 
dimension from east to west in km; Pixel Area = area of pixel in km 2 ; R rs (670) and R rs (645) are applied here as a relative measure of suspended particle load; * = significant 
at adjust p-value; A = significant at 0.05 level; ns = not significant. 
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We evaluated a second MLR algorithm for the MODIS S275.295 product 
(MLR_Fichot), recently published by Fichot et al. (2013), that performed 
slightly better than our MLR algorithm (Table 7). 

3.4. Satellite validation of R rs products 

In order to understand the potential causes of the poor validation 
performance of the QAA-based algorithms, we examined how well the 
satellite R rs compared with the in situ R rs . The validation metrics applied 
to the CDOM properties were also applied to the SeaWiFS and MODIS R rs 
evaluation. The results reveal that the satellite R rs products did not 
match the in situ R rs well in most bands with MAPD ranging from 
29.7% ± 23.1 for SeaWiFS R rs (555) to 102% ± 87 for SeaWiFS R rs (412) 
and from 14.5% ± 10.4 for MODIS R rs (531) to 41.3% ± 26.5 for MODIS 
R rs (412) (Table 8). The band ratios yielded similar results except for 
those based on R rs (490) or R rs (488), which show relatively good agree- 
ment between satellite and field measurements. The other metrics gen- 
erally show a similar pattern with the green satellite R rs bands agreeing 
more closely to the in situ data than the blue and red bands (Table 8). 
One notable exception was in the MODIS average percent bias for 
R rs (412), R rs (443) and R rs (469), which were significantly less than for 
the other bands ( — 2.3 to 3% compared to — 13.8% to — 26.4%). Some 
of these results may be related to the relatively small sample size of 9 
for MODIS-Aqua to 28 matchups for SeaWiFS. Nevertheless, Goyens 
et al. (2013) examined 364 coastal matchups of MODIS-Aqua and in 
situ data and obtained similar results with MAPD of 11-13% for water 
leaving radiances at 488, 531 and 547 nm and 21-35% at 412, 443 and 
667 nm for the same atmospheric correction processing applied in this 
study. 

3.5. Evaluation of factors that could influence satellite product uncertainty 

To explore the possible factors that influence SeaWiFS and MODIS 
retrievals for our study region, Spearman rank correlation analyses 
were conducted on the residuals for satellite-in situ R rs matchups (C sat 
— Q n situ) with the many parameters that could influence the 
uncertainty of satellite R rs retrievals such as solar or sensor geometry, 
atmospheric aerosol and trace gas properties, spatial and temporal 
mismatch between satellite overpass and in situ measurements, and 
surface water turbidity (from R rs (670) and R rs (645), for SeaWiFS and 
MODIS, respectively). A statistically significant positive correlation 
would suggest that the satellite R rs is lower than the in situ measure- 
ments at low values of the correlative parameter and greater than in 
situ at higher values of the correlative parameter. For a negative corre- 
lation result, satellite R rs are higher than the in situ measurements at 
low values of the correlative parameter and lower than in situ at higher 
values of the correlative parameter. A significant correlation between 
the residuals and a particular parameter does not constitute a causation 
of the errors in the matchups, but rather at topic to investigate further. A 
subset of these comparisons reveals that atmospheric properties as well 
as solar zenith angle (SZA), sensor zenith angle and pixel resolution are 
highly correlated with the residuals of R rs from SeaWiFS (Table 9). 
Specifically, the residuals for the SeaWiFS matchups are significantly 
correlated to SZA (negative correlation) and to aerosol optical thickness 
(AOT; positive correlation) for all the R rs bands, but only a subset of the 
bands for sensor zenith angle (positively except for ^5(490)^(555)), 
tropospheric column N0 2 concentrations (N0 2 tropo; negative correla- 
tion), and pixel size x (east to west dimension) and pixel area (positive 
correlation). For MODIS, only the tropospheric column N0 2 concentra- 
tions (N0 2 tropo) and epsilon are significantly correlated to the R rs 
residuals, with the exception of SZA for R rs (678) and time difference 
for R re (443) (Table 9). The small sample size likely limits the robustness 
of the MODIS R rs correlation analysis. 

The Spearman rank correlation analysis was also conducted on the 
residuals of the validation matchups from each of the CDOM and spec- 
tral slope algorithms. For the SeaWiFS a CDOM { 355) and Ocdom(412) 
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Ozone = atmospheric total column ozone; N0 2 strat = stratospheric column N0 2 concentration. 


Table 11 

Spearman rank correlation values of MODIS-Aqua Ocdom^ 12, o C dom355, S 275.295 , and S 30 o : 6oo validation residuals versus various MODIS parameters (atmospheric, sensor geometry, pixel size, etc.). P-value adjusted for number of comparisons, 
p = 0.05/15 = 0.00333. See Table 9 for further details. 
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algorithms, the residuals are negatively correlated with N0 2 tropo, A0T 
and R rs (670), but positively correlated with the angstrom exponent 
(Table 10). The results are similar for all of the a CD0M algorithms evalu- 
ated. The correlations for the SeaWiFS spectral slope algorithms only 
show significant correlations for the MLR S 30 o : 6oo with AOT and the 
QAA-based S275.295 algorithm with N0 2 tropo, AOT and angstrom. 

The correlation analyses for the MODIS CD0M products yield some- 
what different results. The residuals of the QAA-based Qcdom( 355) and 
Qcdom(412) algorithms are negatively correlated with epsilon and 
angstrom whereas the band-ratio and MLR residuals are significantly cor- 
related with SZA (positive correlation) and negatively correlated with 
stratospheric column N0 2 (N0 2 strat), AOT and angstrom (Table 11). 
The correlations with epsilon are also significant and negative for the 
band-ratio algorithm residuals. The QAA-based S 27 s :2 95 residuals are 
only correlated with epsilon at the adjusted P-value threshold, and the 
MLR S 2 75:295 residuals are correlated to SZA (negative), ozone (positive), 
N0 2 strat (positive), AOT (positive), angstrom (positive), and R rs (645) 
(negative). However, neither of the S 30 o:6oo algorithm residuals are signif- 
icantly correlated (at the adjusted p-value) with any of the parameters 
investigated. 

3.5. Extrapolation ofCDOM properties in the UV 
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Information on CD0M in the UV can be extrapolated from visible 
band reflectances because of the strong correlation between o C domM 
with a CD0M at other wavelengths such as a CDOM { 412) (Fig. 8). The 
Pearson correlation analysis between cicdom(412) and c!cdom(250:600) 
of over 1 900 discrete samples from our study region shows that the cor- 
relation coefficient (R) varies from a minimum of 0.76 at 599 nm, 0.89 at 
250 nm and >0.950 between 296 nm and 511 nm (Fig. 8). If we limit our 
data geographically to only the southern MAB and by season (combined 
fall, winter, spring), then the R-value is significantly higher at the lower 
and higher wavelengths (0.99 at 250 nm and 0.97 at 530 nm). 


3.6. Evaluation of satellite-derived CDOM properties 

To further evaluate the applicability of the acDoiw(A). S 30 o: 600 , and 
S 2 75:295 MLR algorithms, temporal and spatial trends of SeaWiFS- 
and MODIS-derived CDOM products were compared with field 
measurements. Several time series of daily MODIS-Aqua and SeaWiFS 
Qcdom(412), S 275:2g 5 , and S 300:600 were generated for four locations 
from the lower Chesapeake Bay, bay mouth, plume and adjacent shelf 
waters that were sampled multiple times from 2004 to 2007. Although 
the sampling dates are not coincident with the satellite observations 



Wavelength (nm) 

Fig . 8 . Pearson correlation coefficient values (R) for the comparison of o C dom( 412) versus 
Qcdom(^) from 250 to 600 nm. Data for all regions represent a compilation of 1933 discrete 
samples collected at various depths within estuaries and continental margin of the 
northeastern U.S. see Fig. lc and Table 1 for geographic distributions and dates of sample 
collections. 
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(with a couple of exceptions), the results show that satellite-derived 
< 3 cdom( 412 ), S 2 75:295, and S 3 oo:6oo (not shown) do capture the general 
trends and variability of CDOM properties (Fig. 9). Despite the fact 
that no data from the lower Chesapeake Bay site were included in the 
algorithm development due to the absence of radiometry data, our algo- 
rithms were able to capture the dynamic range for this site. In addition, 
multi-year mean monthly composites (2004-2007) of satellite 
Qcdom( 412) and S were computed for each of the sub-regions (Fig. f c) 
for comparison with field measurement averaged across each sub- 
region. The monthly composite plots demonstrate that the spatial and 
temporal trends in the satellite data are equivalent to those from the 
field data to within the uncertainty of the validation statistics (Fig. 10; 


only MODIS-Aqua results shown). Note that the sample size for the 
mean monthly field measurements is quite low (Table 2; Fig. 10). 

SeaWiFS and MODIS-Aqua monthly composites of a CD0M ( 412), 
^ 275:295 and S 30 o: 6 oo were processed for multiple years (2004-2007) with 
the better performing MLR and band-ratio algorithms for cicdom(412) 
and MLR and QAA-based algorithms for S 275:295 and S 30 o:6oo to examine 
the performance of algorithms with respect to the spatial and temporal 
variability of these CDOM properties. The processed monthly satellite 
images (only 2006 presented) demonstrate that these algorithms can 
be applied to generate SeaWiFS and MODIS-Aqua maps of these CDOM 
properties (Figs. 11-13). The satellite distributions of a CD0M ( 412) appear 
quite reasonable based on our field measurements and the processes 


a 






Fig. 9. Multi-year time series comparisons of field measurements and (a) daily MODIS-Aqua o C dom(412), (b) SeaWiFS <1cdom(412), and (c) MODIS-Aqua S 2 75-295 for sites located in the 
lower Chesapeake Bay, mouth of the bay, plume, and inner shelf (see Fig. lc). The error bars on the satellite data represent the standard deviation of the 3 x 3 satellite pixel array centered 
on the field station location. Note that field measurements do not coincide with the day and time of the satellite data (with a few rare exceptions). On a few occasions, multiple samples 
were collected at a station over the course of 45 min to 1 h, which is why some plots show three field data points for the same date. 
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Fig. 9 (continued). 


that impact CDOM levels along the continental shelf of the northeastern 
U.S. (Table 2; see Section 4.3). Both field measurements and satellite 
data show a decrease in a CDOM ( 412) along a gradient from nearshore 
to offshore and from north-to-south between the Gulf of Maine and 
southern MAB, excluding Georges Bank (Figs. 10-11; Table 2). 
Qcdom{ 412) is generally higher in spring and fall and lowest during sum- 
mer. The satellite and field data for both S 2 7 s : 2 9 s and S 300 : 6 oo demonstrate 
similar spatial and temporal trends with both S 275.295 and S 30 o: 6 oo decreas- 
ing from nearshore to offshore and from north-to-south (excluding 
Georges Bank) and generally higher values in summer than other seasons 
(Figs. 10b, 12-13; Table 2). For both S 275.295 and S 300 : 6 oo< the SeaWiFS and 
MOD1S QAA-based values are much lower than the MLR values across all 
seasons (Figs. 12-13). 


4. Discussion 

4.1. Algorithm perfoimance 

Using in situ radiometry and CDOM absorption and spectral slope 
measurements, we developed and tuned several types of algorithms 
for OcdomM' S 2 75; 2 95 and S 300 : 6 oo which included several band-ratio 
algorithms, QAA-based, MLR and K d (K) relationships. The in situ data 
show that these algorithms were applicable across the entire study 
region with the exception of the exponential decay band ratio of 
R rs (490)/R rs (555), which required different tuning coefficients for the 
three sub-regions, western Gulf of Maine, New York Bight (region be- 
tween northern New Jersey and Long Island), and southern MAB, 
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Fig. 9 ( continued ). 


where we had both in situ radiometry and CDOM measurements 
(Figs, la and 2a). Since the band ratio algorithms containing the 
412 nm band yielded good relationships across all regions, this would 
suggest that contribution of phytoplankton to R rs (490) somehow differs 
among the three sub-regions yielding the three different band ratio re- 
lationships presented in Fig. 2a. The concentration and composition of 
phytoplankton carotenoid pigments (and potentially chlorophylls b 
and c) could impact R rs (490) in areas where phytoplankton are present 
in appreciable abundances as are typical for the study region. Differ- 
ences in phytoplankton species composition between these regions 
and potentially pigment packaging among the different species are 
the likely causes of the divergent relationships of R rs (490)/R rs (555) 


with o CDO m( 412) (Fig. 2a). Thus, the R rs (490)/R rs (555) band ratio is 
not recommended for retrieval of Qcdom(^)- 

Our results demonstrate several robust in situ bio-optical algorithms 
to quantify OcdomW, $ 275:295 and S 30 o:6oo (Tables 3-6) as well as the 
validation of these algorithms for the SeaWiFS and MODIS-Aqua 
sensors (Table 7; Figs. 6-7). Despite the robustness of the bio-optical 
algorithms illustrated for the field measurements, the satellite retrievals 
of a CD0M (\) for some of the algorithms were not as good as one would 
expect based on our algorithm development results (Tables 3-7). In 
particular, the QAA-based algorithm is one algorithm that seemed to 
fit the field measurements quite well (Fig. 3a; Table 4), but performed 
poorly in the satellite retrievals of CDOM properties (Fig. 7; Table 7). 
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The validation results for the band ratio and MLR algorithms demon- 
strate that these algorithms are appropriate for SeaWiFS and MODIS re- 
trieval of a cdomW and the MLR algorithms for retrieval of S 275:2g5 and 
Sjoo:6oo (Table 7). The QAA-based CDOM spectral slopes were insensitive 
to changes in [a t (412) — a t (443)] at lower S 27 5:295 and S 30 o : boo values. 
For SeaWiFS, the MLR and Exp_R rs (412)/R rs (670) algorithms performed 


better for a CDOM (^) than the other algorithms, and for MODIS, the MLR 
and Exp_R re (412)/R re (547) performed better. Because the band ratio algo- 
rithms lose sensitivity at higher and lower a CD0M (\) concentrations 
(Fig. 2), the MLR algorithm is recommended because it can be applied 
to regions with high and low OcdomM such as estuaries and ocean 
waters beyond the continental margin. The satellite distributions of 
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Fig. 10. Comparison of multi-year (2004-2007) averaged monthly MODIS-Aqua composites (left panels) and field observations (2004-2013; right panels) of (a) cicdom(412) and 
(b) S 2 75:295 by sub-region. The error bars shown on the MODIS data plots represent the standard deviation of the satellite pixels within each sub-region. Sample size for averaged monthly 
field observations was >5. 
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Fig. 10 (continued). 


the QAA-based and MLR algorithms for S 27 5.295 and S 300:600 compared to 
field data indicate that the S 275:295 and S 300:600 MLR algorithms yield bet- 
ter results for these products (Table 2; Figs. 12-13). 

The time series analysis from the lower Chesapeake Bay region 
and monthly composites show that MODIS and SeaWiFS retrievals 
of cicdom{ 412) and S 27 5 : 295 using the MLR algorithms compare to 
within the measured uncertainty of the validation statistics 
(Figs. 9-10; Table 7). Nevertheless, at a few time points and locations 
from the lower Chesapeake Bay region times series, Ocdom(412) field 
measurements were noticeably higher (and S 275:2g5 values lower) 


than the satellite time series trend such as on July 4 (inner shelf), 
September 6 (bay mouth and plume) and November 28 (bay mouth 
and plume) of 2006 and April 23, 2007 (inner shelf). On all but the 
April 23 date, the satellite observations and field sampling occurred 
from one to many days apart. River discharge into Chesapeake Bay 
and freshwater outflow at the bay mouth were significantly higher 
than baseflow in the days and weeks prior to July 4, September 6 and 
November 28, 2006 (see Fig. 11 in Mannino et al., 2008). This could ex- 
plain the much higher in situ Qcdom( 412) (and lower S 275:2g5 ) on these 
dates compared to other sampling periods. The inner shelf sample 
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Fig. 11. Monthly mean satellite-derived a CDOM (412) from SeaWiFS and MODIS-Aqua for the year 2006 computed with the MLR and band ratio (exponential one-phase decay R„(412)/ 
R rs ( 670) orR rs (412)/R„(547)) algorithms. 


from April 23, 2007 was collected -1.75 h and -2.5 h prior to the 
SeaWiFS and MODIS-Aqua overpasses, respectively, but at the bay 
mouth site only -0.5 h after the SeaWiFS overpass and -1 min prior 
to the MODIS overpass. The sensor zenith angle for SeaWiFS (-56°) 
was quite high for the lower Chesapeake Bay region on April 23 
resulting in pixel areas of -4.2 km 2 at the inner shelf site and 
-4.4 km 2 at the bay mouth station, which would likely increase the un- 
certainty in SeaWiFS retrievals. The measured o C dom{ 412) at the bay 
mouth site on April 23 was only 20% higher (7.5% higher for S 275 . 295 ) 
than the MODIS value, which is within the estimated validation uncer- 
tainty (Table 7). The sensor zenith angle (-11°) for MODIS-Aqua and 
pixel areas (-1.3 km 2 ) were quite modest on this particular date. 


4.2. Factors contributing to uncertainty in SeaWiFS and MODIS CDOM and 
R rs products 

SeaWiFS and MODIS R rs products are validated with in situ matchups 
within the NES LME to account for the performance of the CDOM algo- 
rithms. The R rs bands in the blue (412 and 443 nm) and red (667, 670 
and 678 nm) had higher relative error (MAPD) than the other bands 
for both SeaWiFS and MODIS (Table 8). Our R rs matchup results are con- 
sistent with previous studies comparing matchups of SeaWiFS and 
MODIS-Aqua with in situ water-leaving radiances (or R rs ) within marine 
coastal waters (e.g., Antoine et al., 2008; Goyens et al., 2013; Werdell, 
Franz, & Bailey, 2010; Zibordi, Berthon, Melin, D'Alimonte, & Kaitala, 
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Fig. 12. Monthly mean satellite-derived S 300:600 from SeaWiFS and MODIS-Aqua for the year 2006 computed with the MLR and QAA-based algorithms. 


2009). These studies and others affirm that MODIS and SeaWiFS R rs re- 
trievals have higher uncertainties in coastal waters compared to pelagic 
waters of the global ocean (Bailey & Werdell, 2006). 

There are several possible explanations for the discrepancies be- 
tween satellite and in situ R rs measurements and derived products: 
(1) uncertainties in the field-derived measurements (2) in-water con- 
stituents and thus optical properties of coastal waters are too complex 
and variable spatially and temporally for proper satellite-to-in situ 
comparisons (case 2 waters; IOCCG 2000), (3) imperfect atmospheric 
correction of the atmospheric properties associated with the scattering 
and absorption by aerosols and absorption by ozone and nitrogen 
dioxide (N0 2 ) and other atmospheric constituents that can occur in 
higher concentrations on continental margins than the open ocean, 


(4) uncertainties in satellite retrievals due to solar illumination and 
sensor viewing geometry including bidirectional reflectance distribu- 
tion function (BRDF), (5) uncertainties in the satellite sensor calibration 
and vicarious calibration, and (6) mismatch in satellite pixel and sample 
location. All of these factors are likely to impact satellite retrievals of R rs 
and CDOM in our study region to varying degrees. The band ratios do 
yield improved validation results compared to the individual R rs bands, 
particularly for SeaWiFS (Table 8). This suggests some spectrally depen- 
dent uncertainties in the satellite and/or field data R rs bands are 
removed when one band is divided by another band as observed in 
prior studies (e.g„ Zibordi et al„ 2009). 

The Spearman rank correlation analyses of the satellite and in situ 
matchup residuals indicate that certain factors including solar 
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Fig. 13. Monthly mean satellite-derived S 275:2g5 from SeaWiFS and MODIS-Aqua for the year 2006 computed with the MLR and QAA-based algorithms. 


illumination, sensor geometry and atmospheric corrections for aerosols 
and trace gases should be explored further as potential sources of 
uncertainty in R rs and the derived satellite products such as CDOM 
(Tables 9-11 ). The default BRDF corrections within SeaDAS may not 
be appropriate for the case 2 waters within our study region (Morel, 
Antoine, & Gentili, 2002). The SeaDAS default f/Q correction factor 
was applied in the processing of the SeaWiFS and MODIS data used 
for this work because our study region spans both case 1 and case 2 wa- 
ters. The similarity in absorption spectra of CDOM and detrital particles 
may have an impact on our retrievals of CDOM products in more turbid 
areas of our study region (Pan et al., 2008). The Spearman correlation 


results reveal a significant correlation between residuals of q C d om and 
R rs ( 670), a relative measure of particle load, for SeaWiFS but only for 
the band ratio and QAA-based algorithms and not for the MLR algo- 
rithms or the S products (Table 10). For MODIS, only the S 275.295 MLR al- 
gorithms have a significant correlation with particle load, R rs (645) 
(Table 1 1 ). It is noteworthy that the difference in time between the sat- 
ellite overpass and in situ sample collection and satellite pixel dimen- 
sions are not significantly correlated with any of the cicdom(412), 
S 275:295 and S 300 : 6 oo residuals. If spatial variability and temporal variabil- 
ity are major factors in satellite retrieval uncertainties for our CDOM 
products, then this is likely to occur at sub-pixel scales in addition to 
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spatially variable regions such as riverine and estuarine plumes, fronts 
and eddies. Additional validation datasets and satellite data analysis 
would be required to resolve which factors have the greatest effect on 
satellite R rs and CDOM product uncertainties. 

4.3. Spatial and temporal variability ofa C DOM(4!2), S 2 75:295 and S 30 o:6oo 

Satellite-derived distributions of <3 cdom( 412), S275.295 and S 30 o:6oo 
are spatially and temporally consistent with field measurements 
(Figs. 9-13; Table 2). Because terrestrial DOM is the main source of 
CDOM to estuaries and coastal ocean, a CD0M ( 412) values should gener- 
ally decrease from the estuary to the edge of the continental slope, espe- 
cially during the high river flow periods of January to April (Del Vecchio 
& Blough, 2004; Mannino et al., 2008) as illustrated in our satellite im- 
agery and field data (Figs. 10a and 11; Table 2). The MLR a CDOM (412) 
products are more consistent with our expectations than the band- 
ratio algorithms, particularly for inner shelf and estuarine regions 
where o C dom( 412 ) should be higher than mid- and outer-shelf regions. 
In contrast, the band-ratio values are very similar between the lower 
estuaries and inner- and mid-shelf during the winter-spring period. 
The decrease in a CDOM (412) across the shelf from April through August 
is consistent with solar photo-bleaching of CDOM for this region (Del 
Vecchio et al., 2009). The increase in surface a CDO M(412) from early to 
late autumn corresponds to greater river discharge compared to 
summer and seasonal mixing of the water column, which introduces 
non-photo-bleached CDOM into the surface waters from below the 
mixed layer. SeaWiFS and MOD1S distributions of S 27 5:295 and S 300 . 6 oo 
are consistent with the general interpretations of o C dom( 412) distribu- 
tions that seasonal river discharge of terrestrial DOM, summer photo- 
bleaching and seasonal stratification and mixing of the water column 
regulate the seasonal variability of CDOM along the northeastern U.S. 
(Figs. 10-13). In addition to seasonal changes, o C dom( 412), S275.295 and 
S 30 o: 60 o fluctuate quite rapidly on the scale of < 1 h to a few days in coast- 
al waters implying that CDOM could be used to follow dynamic physical 
processes such as freshwater discharge, coastal currents, upwelling, 
eddies, storm events, etc. (Fig. 9). 

4.4. Potential of UV radiometry bands for satellite retrieval of CDOM 
properties 

Satellite observations from SeaWiFS and MODIS are limited to re- 
trievals of ocean reflectance in the visible and near-infrared wavelength 
spectrum. Future NASA ocean color satellite sensors such as on the PACE 
and GEO-CAPE missions are expected to be equipped with UV capability 
and designed to produce improved satellite retrievals for blue bands 
such as 412 nm (Pre-aerosol, clouds, and ocean ecosystem (PACE) 
mission science definition team report, 2012). Other planned sensors 
such as the Japanese Second Generation Global Imager (SGL1) and Geo- 
stationary Ocean Color Imager II (GOC1-II) have planned capabilities to 
retrieve ocean reflectances in the UV. Since CDOM absorbance is more 
intense further into the UV in combination with a weaker absorbance, 
with respect to CDOM, from phytoplankton (pigments) and detrital/ 
non-pigmented particle absorbance from 412 to 350 nm, the addition 
of UV bands should improve the retrievals of CDOM absorption and sep- 
aration of CDOM from phytoplankton and detrital particle absorption. 
For example, prior work has demonstrated the potential for more robust 
algorithms to retrieve CDOM or CDM with UV R rs from a band ratio of 
325/565 and 340/565 (Tedetti et al., 2010) or with UV I< d from a band 
ratio of 320/780 (Hooker, Morrow, & Matsuoka, 2013). Our application 
of UV R rs bands shows promise for improved retrievals of a CD0M , S 275.295 
and S 30 o:boo using the UV-MLR with R rs (380) or K d -based algorithms 
(Tables 5-6). These bio-optical algorithms yield a much better fit to 
the in situ data with R rs (380) and I( d ( 340) or K d (380) through the full 
range of our a CD0M , S 2 7 5:295 and S 30 o:6oo measurements (Figs. 4-5). In ad- 
dition, UV bands may improve the capability of models that retrieve 
multiple inherent optical properties (IOP) (Lee et al., 2002; Maritorena 


et al., 2002; Tilstone et al., 2012; Werdell et al., 2013) to distinguish 
between CDOM and detrital particle absorption. 

4.5. Implications 

Because CDOM is a dominant light-absorbing constituent in the 
ocean and has absorption spectra that overlap significantly with phyto- 
plankton, detritus and minerals, the impact of CDOM absorption on 
water-leaving radiances should be considered when applying ocean 
color remote sensing data to quantity other ocean constituents. Indeed, 
the overlap in absorption results in substantial uncertainties in satellite 
retrievals of Chi (e.g., Pan, Mannino, Russ, Hooker, & Harding, 2010; 
Sauer, Roesler, Werdell, & Barnard, 2012; Siegel, Maritorena, Nelson, & 
Behrenfeld, 2005; Siegel, Maritorena, Nelson, Behrenfeld, & McClain, 
2005; Siegel et al., 2013) using the standard band ratio algorithms 
(OC4 and OC3 for SeaWiFS and MODIS, respectively; O'Reilly et al., 
2000). The capability to retrieve a cdomW and S 300 600 relatively accu- 
rately with the MLR algorithms can enable improved retrievals of 
Chi and potentially other constituents. For portions of our study re- 
gion where detrital particle absorption is very low or insignificant, 
which is generally the case for much of the offshore waters of our 
study region (Pan et al., 2008; Yentsch & Phinney, 1997), existing 
semi-analytical models (e.g., GSM01; Maritorena et al., 2002) can be 
modified to incorporate satellite-derived Qcdom(^) and S 30 o : boo to repre- 
sent CDM and permit retrieval of Chi (and phytoplankton absorption) at 
higher accuracy than the standard maximum band ratio algorithms. 

SeaWiFS and MODIS sensors lack the capability to retrieve ocean 
reflectance in the ultraviolet region of the electromagnetic spectrum. 
Yet, prior work (Fichot et al., 2013; Johannessen et al., 2003; 
Mannino et al., 2008; Swan, Nelson, & Siegel, 2013) and results pre- 
sented here demonstrate retrievals of a CDOM or S for UV wavelengths. 
Because of the strong correlation between o C dom(^) with a CDO M at 
other wavelengths such as o CDO m( 412) (Fig. 8), one can infer infor- 
mation on CDOM in the UV from visible band reflectance. Therefore, 
we can apply SeaWiFS, MODIS and other ocean color sensors with a 
similar band set to retrieve o C dom and S at UV wavelengths. One advan- 
tage to this is the enormous potential to applying satellite retrievals of 
Qcdom(UV), S 30 o:6oo and S275.295 to retrieving chemical properties of 
DOM (e.g., Heines & Benner, 2003; Mannino et al., 2008; Spencer, 
Aiken, Wickland, Striegl, & Hernes, 2008; Spencer et al., 2010; 
Spencer, Butler, & Aiken, 2012; Fichot & Benner, 2011; Fichot & 
Benner, 2012; Fichot et al., 2013). 

Since the magnitude of a CDOM (^) and CDOM spectral slopes depend 
on its origin and composition, satellite retrievals of CDOM properties in 
coastal waters can be applied to study relative contributions of different 
carbon sources, track water masses, quantify distributions and 
variability of DOC, and estimate photooxidation of DOM. Satellite 
products of Ocdom(R) can be applied to retrieve DOC concentrations 
(e.g., Griffin, Frey, Rogan, & Holmes, 2011; Korosov, Posdnyakov, & 
Grassl, 2012; Mannino et al„ 2008) and quantify processes such as 
river export of DOC (Del Castillo & Miller, 2008; Lopez, Del Castillo, 
Miller, Salisbury, & Wisser, 2012), cross-shelf fluxes of DOC, and photo- 
chemical production of inorganic carbon (Belanger et al., 2008; Del 
Vecchio et al., 2009; Fichot & Miller, 2010). CDOM can be used as a 
good proxy for tracing physical circulation and water-mass history, 
tracking freshwater plumes, identifying pollutants and assessing their 
impact on inland and coastal water quality and ecosystem health. 
S 30 o:6oo can vary with the source of DOM, but S 300:600 is also influenced 
by photochemical and biological processes (see reviews by Blough & 
Green, 1995; Blough & Del Vecchio, 2002). Since S 300:6 oo generally in- 
creases with salinity and decreasing molecular weight and aromatic 
content (Blough & Del Vecchio, 2002), this suggests that lower S 30 o:6oo 
values represent higher proportions of terrestrial DOM, though there 
are notable exceptions with respect to the salinity relationship 
(e.g., Stedmon et al., 2000). Helms et al. (2008) recently introduced 
S 275:295 as a more powerful measure of DOM source, molecular weight 
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and of photochemical and microbial modifications to CDOM than S 300:60 o. 
With validated algorithms for Qcdom(^). $ 275:295 and S 30 o:6oo such as those 
presented here, it is possible to retrieve other chemical properties of DOM 
through in situ relationships between OcdomOs) or S275.295 and dissolved 
lignin phenol concentrations (e.g., Hernes & Benner, 2003; Spencer 
et al„ 2008, 2009, 2010, 2012; Osburn & Stedmon, 2011), lignin yields 
(Fichot & Benner, 2012; Fichot et al„ 2013), and dissolved black carbon 
(e.g., Stubbins, Niggemann, & Dittmar, 2012). 

5. Conclusions 

Several algorithms have been developed and validated for the re- 
trieval of a CDOM (\), S 275:295 and S 300 :6oo with SeaWiFS and MODIS- 
Aqua for the continental margin of the northeastern U.S. The validation 
metrics and processing of SeaWiFS and MOD1S data illustrate that the 
MLR algorithms provide the optimal retrievals for OcdomW, S 30 o:6oo 
and S 27 5 : 295 . The relative errors (MAPD) for the validated CDOM prod- 
ucts fall within the uncertainty requirements of 35% for satellite retriev- 
al of open ocean chlorophyll-a (McClain, 2009). The MAPD for a CDOM (\) 
algorithms range from -24% to 30% and from 6% to 10% for the S275.295 
and S 300 :soo products, respectively. Algorithms developed and evaluated 
with in situ measurements employing R rs bands in the ultraviolet 
produce better fits to the in situ data than the algorithms with only 
visible wavelength bands. The algorithms with UV bands will enable 
improved retrievals of a cdoim(A). S 2 75:295 and S 30 o:6oo from future 
satellite sensors with UV measurement capabilities. The correlation 
analyses of the validation residuals for the CDOM products and R rs 
suggest that several factors should be explored further to understand 
and possibly diminish the uncertainty in satellite retrievals within the 
coastal ocean, including solar zenith angle, sensor viewing angle, 
atmospheric products applied for atmospheric corrections, as well as 


the in situ spatial and temporal heterogeneity of a CD0 M and other optical 
properties. 
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Appendix A 


Table Al 

Coefficients and goodness of fit parameters for band-ratio bio-optical a CDO M algorithms developed from field observations of o CDO m(^) and remote sensing reflectance (R rs (\) ). See Table 2 
for details. 


Exponential one-phase decay equation 

Qcdom 275 

Gcdo/w355 

Qcdom 380 

Qcdom4 12 

Qcdom443 

X = Ln[(Y — B 0 ) / B 2 ] / — B, 

B 0 = 0.2581 

0.2452 

0.2492 

0.2487 

0.2479 

Y = R rs (412)/R rs (555) 

Bt = 1.583 

5.576 

8.689 

14.028 

23.40 

X = Qcdom(X) 

B 2 = 24.87 

4.838 

4.608 

4.085 

3.770 

R 2 

0.881 

0.942 

0.948 

0.950 

0.937 

Sy.x 

0.113 

0.086 

0.081 

0.0806 

0.084 

n 

142 

151 

152 

153 

149 

Minimum band ratio 

0.31 

0.295 

0.295 

0.295 

0.295 

X = Ln[(Y — B 0 ) / B 2 ] / — B, 

B 0 = 0.9925 

0.8569 

0.865 

0.8625 

0.8502 

Y = R„(412)/R rs (667) 

Bi = 2.054 

7.661 

11.55 

18.44 

30.53 

X = Qcdom(X) 

B 2 = 634.2 

91.97 

79.16 

62.89 

54.78 

R 2 

0.845 

0.955 

0.959 

0.943 

0.938 

Sy.x 

1.046 

0.720 

0.694 

0.736 

0.788 

n 

139 

144 

145 

147 

146 

Minimum band ratio 

1.29 

1.1 

1.1 

1.1 

1.1 


Table A2 

Satellite validation statistics of CDOM exponential decay band ratio algorithms for SeaWiFS and MODIS-Aqua. 

Product 

Algorithm 

MAPD 

±sd 

RMSE 

Median ratio 

SIQR 

%Bias 

Slope 

R 2 

n 

Qcdom275 

Exp_R„(412)/R rs (555) 

24.8 

22.3 

1.77 

0.88 

0.209 

-27.3 

0.202 

0.45 

41 

Qcdo/vj355 

Exp_R„(412)/R rs (555) 

32.5 

22.0 

0.348 

0.73 

0.222 

-34.9 

0.282 

0.49 

41 

a CDOM^ 80 

Exp_R„(412)/R rs (555) 

33.3 

21.5 

0.224 

0.68 

0.206 

-36.4 

0.289 

0.50 

41 

Qcdoat4 12 

Exp_R rs (412)/R rs (555) 

35.6 

20.3 

0.130 

0.62 

0.199 

-37.6 

0.313 

0.51 

41 

Gcdom443 

Exp_R„(412)/R rs (555) 

36.6 

20.3 

0.0746 

0.61 

0.205 

-38.1 

0.328 

0.52 

41 

q cdom275 

Exp_R„(412)/R rs (667) 

28.7 

22.8 

1.036 

1.16 

0.176 

6.6 

0.287 

0.24 

46 

Qcdo/vj355 

Exp_R„(412)/R rs (667) 

31.6 

25.0 

0.196 

1.20 

0.219 

4.9 

0.350 

0.35 

46 

Qcdom380 

Exp_R„(412)/R rs (667) 

29.6 

23.0 

0.122 

1.17 

0.216 

3.2 

0.378 

0.37 

46 

Qcdom^ 12 

Exp_R„(412)/R rs (667) 

27.6 

18.2 

0.0683 

1.14 

0.225 

0.19 

0.429 

0.41 

46 

Gcdom443 

Exp_R„(412)/R rs (667) 

27.8 

17.6 

0.0391 

1.13 

0.227 

0.16 

0.448 

0.41 

46 
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